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20.  Abstract  (continued) 

an  average  value  of  2  sec  and  (like  the  previously-measured 
energy  conversion  yields)  appears  to  increase  with  the  height 
of  the  maximum  of  the  energy  deposition  rate  profile.  Fits  of 
radiance  data  from  moving  arcs  to  a  simple  model  of  vibrational- 
bands  emission  from  nitric  oxide  result  in  exponential  predosing 
time  constants  between  1£  and  2 i  sec,  in  qualitative  agreement 
with  current  aerochemi stry  models. 

Volume  emission  rate  distributions  of  the  blue-violet  air 
fluorescence  within  20  m  of  EXCEDE's  accelerator  were  unfolded 
from  radiometric  photographs  taken  at  characteristic  ejection 
altitudes  from  onboard.  This  emission  decreases  very  rapidly 
away  from  the  rocket,  and  transverse  to  the  beam  axis  fits  a 
gaussian  whose  width  first  increases  and  then  decreases.  The 
measured  excitation  patterns,  which  cannot  be  fully  attributed 
to  outgassing  from  the  vehicle,  are  interpreted  in  terms  of  a 
recent  theoretical  treatment  of  discharge  initiation  and 
electron  heating  resulting  from  an  interaction  between  the 
charged  beam  and  the  local  plasma  it  creates.  The  effects  of 
such  a  collective  process  on  feature  radiances  measured  on 
upleg  and  downleg  are  discussed,  and  the  altitude  dependence  of 
distance  along  the  injected  beam  beyond  which  the  secondary 
exciting-electron  component  no  longer  measurably  affects  the 
infrared  radiation  spectrum  is  estimated. 


SUMMARY 


The  objective  of  the  work  reported  here  is  to  inter¬ 
pret  rocket  and  aircraft  data  on  the  infrared-optical  radia¬ 
tions  resulting  from  naturally-occurring  and  artificially- 
injected  energetic  particles  that  simulate  excitation  by 
atmospheric  nuclear  explosions.  Photometrically-calibrated 
images  from  onboard  of  EXCEDE:  Spectral's  ejected  electron 
beams  and  high-resolution  auroral  radiometry  data  from  the 
Air  Force  Geophysics  Laboratory's  instrumented  aircraft  are 
considered.  The  field  programs,  performed  by  AFGL  under 
Defense  Nuclear  Agency  sponsorship,  are  described  in  Ref's 
1,2,5  (EXCEDE)  and  28  (the  Test  Plan  for  the  aircraft-borne 
measurements) . 

The  correlation  between  2.8  -  3.1  urn  infrared  emission 
from  and  particle  energy  input  to  the  E  region  was  systemati¬ 
cally  determined  from  fifteen  of  the  highest  signal/  noise 
segments  of  data  taken  in  or  near  the  magnetic  zenith  in 
>IBC  II  aurora  by  coaligned  radiometers  having  <1  km  instan¬ 
taneous  projected  footprint.  The  results  apply  to  air  parcels 
typically  30  km  in  height  centered  at  altitudes  determined 
by  two-color  photometry  from  the  aircraft.  Application  of  a 
standard  cross-correlation  method  to  the  input  and  output 
data  streams  resulted  in  peaks  in  the  correlation  function 
that  cluster  near  2  sec,  with  individual  values  ranging  between 
0  and  7  sec.  Like  the  previously  determined  (in  Ref  1)  energy 
conversion  yield  under  steady-state  bombardment,  these  lags  in 
the  SWIR  emission  give  some  indication  of  increasing  with  the 
height  of  the  maximum  of  the  profile  of  energy  deposition  by 
auroral  electrons.  This  behavior  and  the  magnitude  of  the 
delay  are  in  qualitative  agreement  with  the  current  model  of 
excitation  of  overtone  vibrational  chemiluminescence  from  NO 
molecules  by  energetic  particles  from  nuclear  explosions. 


Fits  of  a  simplified  input-output  model  were  made  to 
three  data  segments  from  isolated  moving  arcs,  to  whose 
rapidly  varying  energy  input  the  results  of  the  procedure 
are  most  sensitive.  They  indicated  that  the  effect  of  pre¬ 
dosing  on  instantaneous  SWIR  yield  decreases  exponentially 
with  time  constants  1-1/2  to  2-1/2  sec,  but  failed  to  resolve 
the  expected  dependence  of  grow-in/decay  time  on  mean  energy 
deposition  altitude.  The  lags  show  directly  the  temporal 
smearing  of  2.8  -  3.1ym  sky  backgrounds  relative  to  par¬ 
ticle  energy  input  (whose  distribution  at  past  high  altitude 
nuclear  bursts  is  readily  determined  from  photographs)  ,  and 
only  somewhat  less  directly  the  attenuation  of  high  spatial 
frequency  components  of  the  atmosphere's  infrared  radiance 
due  to  its  wind  shear  and  turbulence. 

Volume  emission  rates  of  the  blue-violet  fluorescence 
of  the  air  within  20  m  of  EXCEDE's  electron  accelerator  were 
unfolded  from  the  radiances  measured  from  calibrated  black- 
and-white  photograhs  taken  onboard  the  vehicle,  for  the  pur¬ 
pose  of  assessing  the  effects  on  the  spectrums  of  collective 
interactions  of  the  electron  beam  and/or  outgassing  from  the 
rocket.  78  millisec-exposure  images  from  characteristic 
upleg  and  downleg  altitudes,  at  the  fixed  injection  current 
of  7  -  8  amperes,  were  selected  for  the  evaluation.  It  was 
found  that  the  volume  emission  distributions  transverse  to 
the  beam's  axis  can  be  fit  to  adequate  precision  by  gaussians 
The  pattern  within  ~5  m  of  the  accelerator  is  almost  independ 
ent  of  ambient  air  density  between  128  km  apogee  and  110  km 
(and  perhaps  to  lower  altitudes);  more  radiation  than  pre¬ 
dicted  from  independent-particle  transport  theory  is  present 
to  at  least  20  m,  near  which  the  initially-diverging  beams 
appear  to  be  refocusing  (if  this  interpretation  is  correct 
they  do  not  experience  severe  velocity  spreading);  and  the 
volume  emission  rates  are  higher  on  upleg  than  downleg,  more 


so  with  decreasing  rocket  altitude  (or  increasing  time  between 
the  measurements) . 

The  spatial  and  altitude  distributions  cannot  be  fully 
attributed  to  outgassing  of  N2  molecules  from  EXCEDE's  pay- 
load  and  recovery  system.  We  interpreted  them  following  a 
recent  preliminary  model  (Ref  20)  of  electron  heating  and 
initiation  of  discharges  by  an  interaction  between  the  beam 
and  artificially-generated  plasma.  Intensities  of  ultraviolet, 
visible,  and  infrared  features  measured  on  upleg  and  downleg 
(not  all  of  which  have  been  extracted  from  the  data)  are  shown 
to  provide  a  further  check  on  this  model.  The  distances  along 
the  injection  flux  tube  at  which  the  electrons  created  and 
heated  by  the  close-in  discharge  have  little  effect  on  air's 
infrared  emission  spectrum  are  estimated,  and  the  value  for 
125  km  is  found  to  be  in  reasonable  agreement  with  a  ground- 
based  measurement  of  the  beam's  radiance  distribution.  These 
ranges,  whose  precision  can  be  readily  improved  by  more  re¬ 
fined  calculations  of  the  initial  spectrum  and  transport  of 
this  secondary  excitation  source,  apply  in  designing  EXCEDE- 
concept  simulation  experiments  in  which  the  infrared  radia¬ 
tion  outputs  are  measured  from  remote  (daughter)  rocket 
platforms . 


PREFACE 


This  is  a  report  of  PhotoMetrics '  second  year's  work 
under  Contract  DNA001 -81 -C-0003  on  analysis  and  interpreta¬ 
tion  of  atmospheric  nuclear  ef fects-simulation  data  taken 
for  the  Defense  Nuclear  Agency.  The  data  are  from  the 
EXCEDE:  Spectral  electron-injection  rocket  (A51.970), 
launched  19  Oct  79  at  Poker  Flat  Research  Range,  AK ,  and 
flights  at  auroral  latitudes  of  the  Air  Force  Geophysics 
Laboratory's  NKC-135A  infrared-optical  aircraft  55-33120,  in 
1979-80.  Earlier  results  from  these  field  experiments  are 
reported  in  Ref's  1  (the  report  of  the  first  year's  work  on 
the  program),  2,  and  27. 

The  work  was  done  under  the  direction  of  I.L.  Kofsky. 
The  authors  express  their  thanks  to  R.R.  O'Neil,  E.R.  Huppi 
(AFGL) ,  D.R.  Archer  (Mission  Research  Corp.) ,  J.  Barrett,  G.B 
Parrent  Jr  (PhotoMetrics) ,  and  particularly  K.  Papadopoulos 
(University  of  Maryland)  for  providing  useful  information, 
and  to  Mrs.  C.C.  Rice  for  typing  the  manuscript.  The  en¬ 
couragement  of  A.T.  Stair  Jr  of  AFGL  and  LtCol  W.  McKechney 
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SECTION  I 


VOLUME  EMISSION  RATES  OF  AIR  FLUORESCENCE 
FROM  EXCEDE:  SPECTRAL 

BACKGROUND 

Evaluation  to  date  of  the  onboard  (Ref's  1,2)  and  ground- 
based  (Ref's  3,4)  photographs  of  the  airglows  excited  by  the 
3  kV  electron  beam  from  EXCEDE:  Spectral  (A51.970,  19  Oct  79; 
Ref's  5,6)  has  shown  their  radiance  distributions  to  be 
markedly  different  from  those  calculated  assuming  the  energy 
deposition  is  by  independently-acting  particles  (Ref  7) ,  both 
within  individual  scenes  and  in  dependence  on  the  upper  atmo¬ 
sphere's  density.  This  finding  has  been  interpreted  as  due 
to  electric  discharges  from  collective  interaction  of  the 
charged  sounding  rocket  and/or  cold  electron  stream  with  the 
ambient  and  artif icially-created  ionospheric  plasma  (Ref  1 , 
on  the  basis  of  Ref's  8  and  9).  An  important  result  of  this 
interaction  is  that  the  number  densities  and  energy  distri¬ 
butions  of  the  secondary  electrons  in  the  fields  of  view  of 
infrared  instruments  on  the  vehicle  differ  substantially 
from  those  produced  by  individual  charged  particles  (Ref's 
10,11),  such  as  are  emitted  from  atmospheric  nuclear  explosions 
(e-rays,  debris  ions,  secondary  electrons  of  thermal  x-rays 
and  x-rays).  As  these  secondary  electrons  efficiently  excite 
vibrational  states  of  pre-existing  molecules  (principally, 

CO2  and  its  precursor  N2,  and  H2O)  and  also  play  a  part  in 
the  aerochemical  processes  that  create  new  infrared-active 
species  (NO,  O3 ,  N0+ ,  others),  the  difference  in  energy  spectra 
would  introduce  errors  in  scaling  EXCEDE' s  radiometry  data 
to  infrared  sky  background  intensities  in  nuclear  warfare. 

The  presence  of  enhanced  fluxes  of  secondaries  in  the 
beam  was  further  indicated  by  unexpectedly  high  radiances  in 


several  of  the  ultraviolet,  visible,  and  infrared  molecular 
band  systems  of  air  that  have  high  optical  cross-sections  for 
excitation  by  low-energy  electrons,  as  was  detailed  in  Ref  1. 
(While  some  aspects  of  the  spatial  and  altitude  dependence  of 
photographic  brightness  might  appear  due  to  molecules  outgassed 
from  the  rocket  body  --  on  which  more  presently  --  outgassing 
can  not  explain  the  spectral  ratios  of  N2  emission  features.) 
Furthermore,  qualitatively  similar  vis ible -rad iat ion  patterns 
have  been  reported  from  other  injections  of  ~1 —  1 0  keV 
electrons  into  both  the  thermosphere  and  low-pressure  air  in 
laboratory  tanks,  as  also  discussed  in  Ref  1  (more  recent 
observations  are  in  Ref's  12  and  13).  The  optical  data, 
along  with  data  on  radiofrequency-wave  emission  (oscillating 
electric  fields  heat  the  electrons)  ,  plasma  density  and 
temperature,  and  neutralization  of  spaceborne  accelerators, 
have  brought  forth  a  body  of  thought  on  collective  interac¬ 
tions  of  the  charged  beam  and  rocket  with  the  surrounding 
atmospheric  plasma  (as  reviewed  in  Ref's  9,  13,  and  14;  see 
also  Ref  15). 

The  purpose  of  analyzing  EXCEDE:  Spectral’s  radiometric 
photographs  can  be  stated  as  determining  the  distance  from 
the  ejection  rocket  beyond  which  the  effect  of  secondary 
electrons  from  any  collective  processes  on  the  infrared  chemi¬ 
luminescence  and  fluorescence  spectrums  can  be  neglected  or 
reliably  corrected  for.  That  is,  the  work  reported  here  is 
intended  to  characterize  the  plasma  heating  (if  any  occurs) 
only  to  the  extent  required  to  assess  its  interference  with 
the  simulation  of  nuclear  particle-excited  sky  background 
radiations  and  processes.  This  information  about  spatial 
distributions  of  short-wavelength  visible  air  fluoresence  is 
intended  1)  for  interpreting  existing  EXCEDE:  Spectral  optical/ 
IR  data  and  2)  for  designing  experiments  in  which  the  infrared 
sensors  view  excitation  of  particle  beam-excited  air  from 


daughter  platforms  spatially  separated  from  the  mother 
(accelerator-carrying)  vehicle. 

Up  to  now,  evaluation  of  these  photographic  images  (Ref's 
1,2,13  (p217  ff))  has  consisted  of  determining  the  glow's 
surface  radiance  distribution  as  projected  to  the  wide-angle 
camera,  from  five  rocket  altitudes.  In  the  present  work,  we 
have  systematically  computer-reduced  images  at  a  series  of 
upleg  and  downleg  altitudes  with  a  fixed  injected  current  to 
volume  emission  rates  within  -20  meters  from  the  electron 
accelerator.  These  represent  the  "fundamental"  physical 
quantity  in  the  film  data,  as  they  measure  the  volume  rates 
of  excitation  of  N2  molecules  by  beam  and  secondary  electrons 
and  thus  the  sum  of  products  of  excitation  cross-section  and 
flux  of  these  particles.  The  longitudinal  variation  of  this 
excitation  measures  the  extent  of  the  interaction-and-discharge 
region,  and  its  width  transverse  to  the  beam  axis  provides 
further  information  on  the  behavior  of  the  beam  and  plasma. 

INSTRUMENTATION  AND  QUALITATIVE  DESCRIPTION  OF  THE  DATA 

The  rocket's  photographic  system  of  one  16  mm  camera 
filtered  for  sensitivity  to  blue  and  violet  fluorescent  radia¬ 
tion  coaligned  with  a  similar  unit  with  3-layer  color  film  (in 
effect,  Ektachrome)  is  described  in  Ref  1;  for  completeness, 
we  review  the  fast  monochrome  camera's  pertinent  character¬ 
istics  and  calibration  here.  Figure  1  shows  its  spectral 
response.  Figure  2  the  diffuse  transmission  density-log 
exposure  (H&D)  characteristic  of  its  black-and-white  film 
measured  with  the  same  short-pass  filter  (Schott  BG-12  colored 
glass)  in  the  sensitometer  as  was  used  in  the  camera,  and 
Figure  3  the  intensity  vignetting  and  geometric  distortion 
of  its  lens.  The  camera  images  a  62°  x  94°  field  onto  the 
7.64  x  12.0  mm  film  frame,  with  the  long  dimension  in  the 
direction  of  the  ejected  electron  beam  (refer  to  the  photo- 
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Figure  1. 


Spectral  response  of  the  EXCEDE  camera  and 
film.  The  sharp  falloff  above  4600  A  results 
from  the  attenuation  by  its  BG-12  Schott 
filter,  and  its  response  at  wavelengths 
below  3900  A  is  limited  by  lens  absorption. 
The  lengths  of  the  dotted  lines  are  roughly 
proportional  to  the  relative  intensities  in 
an  EXCEDE  spectrum  at  74  km  rocket  altitude 
(see  text). 


FILM  SENSITIVITY 


MAXIMUM  DENSITY 


measured  in  blue  light. 


DISTANCE  FROM  FRAME  CENTER,  mm 

Measured  relative  irradiance  in  the  image 
(solid  line)  and  field-of-view  half-angle 
(dotted  line)  as  a  function  of  displace¬ 
ment  of  image  points  from  the  field's 
center.  These  lens  vignetting  and  dis¬ 
tortion  data  are  applied  in  converting 
the  film  density  distributions  to  scene 
radiance  distributions. 


graphs  and  isophote  plots  in  Ref's  1  or  13)  and  the  diagonal 
sharply  truncated  beyond  ~105°.  Points  1/2  m  away  from  the 
lens  are  in  focus  to  within  the  image's  ~1/3°  circle  of  con¬ 
fusion.  Continuous  sequences  of  8,  78,  and  800  millisec 
exposures  were  made,  uncoordinated  with  the  programming  of 
EXCEDE's  electron  accelerator  pulses. 

Fig  1  also  shows  the  relative  intensities  of  N2+  First 
Negative  (B*X)  and  N2  Second  Positive  (C+X)  band  inten¬ 
sities  measured  by  EXCEDE's  visible-sensitive  spectrometer  at 
73  km  downleg  altitude  (p  28  of  Ref  5).  Within  the  camera's 
passband  this  published  spectrum  is  sensibly  that  of  intense 
aurora  (Ref  16)  and  of  ionizing  radiations  from  nuclear  ex¬ 
plosions  depositing  energy  at  altitudes  between  55  and  125  km. 
Under  conditions  that  the  glow  excited  by  the  beam  has  this 
emission  spectrum,  about  62%  of  the  actinic  exposure  is  due 
to  the  N2+  3914  S  (0,0)  band,  another  14%  is  from  the  4278  & 
(0,1)  band,  and  16%  is  from  the  N2  Second  Positive  bands. 

Where  the  population  of  electrons  in  the  camera  field  with 
energy  above  the  11.03  eV  threshold  for  exciting  N2's  C  state 
is  enhanced  (relative  to  that  resulting  from  impact  of  primary 
electrons)  ,  the  Second  Positive  system  makes  a  larger  relative 
contribution.  (The  threshold  for  exciting  the  B  state  of  N2+ 
is  18.74  eV.) 

Such  a  change  in  EXCEDE's  emission  spectrum  is  evidenced  by 
the  magenta  color  outward  from  the  apparent  beam  axis,  which 
is  due  to  excitation  of  N2 ' s  violet-blue  C+X  and  red  B+X 
(First  Positive)  triplet-state  bands  by  laterally-scattered 
secondary  electrons.  (Refer  to  the  schematized  cross-sections 
in  Fig  17  of  Ref  2.)  This  emission  of  course  also  overlies 
the  central  whitish-blue  "core”  region,  as  the  beam  Coulomb- 
scatters  secondaries  at  all  azimuths.  We  have  not  yet  analyzed 
the  color  film  (which  has  virtually  the  same  high  signal/ 
noise  and  overall  data  quality  as  the  black-and-white  film) 


to  ascertain  what  information  about  the  secondary-electron 
spectrum  is  contained  in  this  purplish  component  of  the 
photographic  radiance  pattern. 

In  any  case,  since  the  camera  has  a  relatively  wave¬ 
length-insensitive  response  (Fig  1)  only  a  small  change 
in  its  radiometric  calibration  results  from  any  such  changes 
in  the  glow's  emission  spectrum.  As  was  shown  in  Ref  1, 
application  of  the  fluorescence  spectrum  of  air  in  the  E-region 
(which  as  noted  is  closely  similar  to  the  published  low-altitude 
EXCEDE  spectrum)  shows  that  the  effective  exposing  wavelength 
is  about  4000  A  in  the  "core,"  and  the  known  cross-sections 
for  excitation  of  Second  Positive  bands  by  the  secondaries 
indicate  that  in  those  regions  where  the  color  film  is  magenta 
it  is  about  3900  8.  The  camera's  spectral  sensitivity  differs 
by  only  12%  over  this  narrow  wavelength  range.  In  assigning 
the  absolute  calibrations  in  Fig  2  we  corrected  the  film's 
response  to  the  sensitometer ' s  illumination  spectrum  to  its 
response  to  4000  A  photons.  Since  the  onboard  spectro¬ 
meter’s  field  of  view  encompasses  the  core  its  spectrums  give 
the  actual  relative  contribution  of  First  Negative  and  Second 
Positive  band  emission  to  the  photographic  signal  at  all 
ejection  altitudes,  in  that  central  area  of  film  image  where 
the  color  remains  substantially  constant.  (These  spectrums 
have  not  yet  been  published.)  In  summary:  the  measured 
total  number  of  photons  emitted  per  unit  volume  of  glow  will 
change  only  slightly  with  its  emission  spectrum,  but  where  the 
fluxes  of  11-  ~19  eV  secondary  electrons  are  enhanced  the 
N2  Second  Positive  system  contributes  more  than  the  afore¬ 
mentioned  fraction  of  the  radiation  signal. 

The  camera's  absolute  radiometric  calibration  is  subject 
to  a  systematic  error  at  low  exposures  due  to  uncertainty  in 
transmission  of  blue  light  by  the  more  optically  dense  area  of 
the  sensitometer ' s  step  tablet.  This  error  is  a  maximum  of  25% 


near  the  toe  of  the  film's  response  characteristic  and  becomes 
negligible  when  the  film-plane  exposure  reaches  1  x  10~1  erg/ 
cm^.  Note  (Fig  2)  that  the  film's  contrast  y  is  essen¬ 
tially  constant  between  the  maximum  exposure  recorded  and  ~10“2 
erg/cm^.  This  latter  exposure  results  from  low  scene  radi¬ 
ances,  of  the  order  of  1/2  erg/cm^-sec-sterad  (1  megarayle igh)  , 
and  as  previous  isophote  plots  (in  Ref's  1  and  2)  show  the 
glow  brightnesses  drop  to  this  level  only  well  outside  the 
flux  tube  into  which  the  electron  beam  is  injected.  Thus 
the • fractional  uncertainty  in  radiance  introduced  by  micro¬ 
densitometry  errors  and  digitization  noise,  which  is  loge10  x 
(density  error)/y,  changes  little  over  most  of  the  glow's 
image.  Where  the  density  error  is  0.05  units,  this  radiometry 
error  would  be  17%. 

A  further  systematic  error  at  all  exposures  is  due  to  un¬ 
certainty  in  the  lens's  transmission  of  axial  rays,  and  a  small 
and  scene  angle -dependent  error  results  from  uncertainty  in  the 
transmission  of  the  blue  glass  filter.  Still  further  radio- 
metric  error  is  introduced  by  imprecision  of  the  correction 
for  vignetting  by  the  camera's  wide  angle  lens  (Fig  3),  which 
becomes  large  toward  the  edge  of  its  field  of  view.  Con¬ 
sidering  these  and  other  potential  errors  noted  in  Ref  1, 
we  estimate  the  absolute  brightness  measurement  accuracy  to 
be  within  a  factor  ±1.5  at  the  geometric  center  of  the 
image,  and  less  than  a  factor  2  at  the  extremes  of  the  frame 
within  ~2  m  from  and  beyond  ~15  from  the  rocket  body. 

Absence  of  any  unforeseen  large  error  in  the  film  radiometry 
is  evidenced  by  the  agreement'  between  the  photographic  radi¬ 
ances  4-1/2  m  down  the  beam  and  those  measured  in  essentially 
the  same  projection  by  a  3914  R-band  photometer  in  the 
payload  (Fig  16  of  Ref  2,  from  yet-unpublished  data). 

Lens  distortion  causes  radiation  originating  at  large 
angles  o  to  image  at  angles  <0  (both  measured  from  the 
principal  points  of  the  lens).  This  small  effect  was  taken 
into  account  by  incorporating  into  the  computer  reduction 


the  measured  relation  (in  Fig  3)  between  position  on  the  film 
image  and  direction  in  object  space.  Offset  of  the  lens  axis 
from  the  geometric  center  of  the  film  frame  can  be  seen  to  be 
less  than  1 50um  in  the  long  dimension  of  the  frame  --  that 
is,  the  direction  of  the  beam  as  is  shown  by  the  position 
of  the  edges  of  the  circular  image.  Misalignment  transverse  to 
the  beam  would  have  only  second-order  effects  on  the  results  of 
the  analysis,  and  insofar  as  the  film  advance  is  pin-registered 
any  photogramme  try  errors  would  be  systematic  throughout  the 
data  set.  Timing  error  is  negligible,  as  each  photographic 
image  can  be  directly  identified  with  accelerator  operation 
conditions  and  rocket  altitude. 

The  camera  viewed  toward  the  electron  beam  from  3-1/4  m 
north  of  the  anode  of  the  accelerator  that  operated  over  most 
of  the  flight,  with  its  optic  axis  in  the  plane  defined  by  the 
injection  field  line  and  the  rocket's  long  axis.  The  projec¬ 
tion  direction  moved  to  21°  east  of  the  meridian  when  the 


rocket's  azimuth  was  changed  just  after  apogee  to  align  the 
payload  with  the  trajectory  so  as  to  increase  its  instruments' 
sight  paths  through  irradiated  air;  this  causes  the  image  of 
the  glow  to  tilt  to  the  right  (northwest)  on  the  downleg 
frames.  A  side  view  of  the  camera's  field,  the  ejected  beam, 
and  the  location  of  these  instruments  --  whose  fields  inter¬ 
cept  the  beam  axis  4-6  m  outward  from  the  accelerator  --  is 
in  Ref  1 ,  and  a  simplified  diagram  appears  in  Figure  4. 

Image  positions  of  some  points  on  the  beam's  axis  are 
1.5  m  from  the  accelerator,  near  the  bottom  edge 
of  the  frame  (limit  of  camera  fov) ; 

1.8  m,  the  line  from  the  camera  lens  is  per¬ 
pendicular  to  the  axis; 

4.4  m,  center  of  the  image,  the  vector  from  the 
lens  makes  a  45°  angle  with  the  axis; 

20  m,  4-1/2  mm  above  the  center  of  the  image, 
the  vector  from  the  lens  makes  a  8-1/2° 
angle  with  the  axis. 


The  direction  at  which  the  finite-width  beam  appears  to  con¬ 
verge  at  infinity  lies  within  the  frame.  We  have  so  far  ex¬ 
tended  the  film  data  analysis  to  20  m  along  the  beam  from  the 
accelerator,  as  the  physical  properties  of  the  glow  do  not  seem 
to  change  substantially  beyond  about  8m--  although  there  is 
now  evidence  of  a  refocusing  node  --  and  the  foreshortened  per¬ 
spective  (and  secondarily,  imprecise  vignetting  correction) 
introduces  increasing  error  into  the  unfolds  of  the  radiance 
distributions.  Specifically,  the  sight  path  that  intercepts 
the  beam  axis  at  20  m  passes  through  the  "edges"  of  the 
emitting  region  at  distances  -10  and  30  m  from  the  accel¬ 
erator;  although  the  column  emission  is  weighted  toward  the 
axis,  this  oblique  projection  may  be  the  reason  that  no 
obvious  beam  refocusing  node  is  seen  in  the  surface  radiance. 
(The  first  node  is  at  21  m  in  the  magnetic  field  over  Poker 
Flat  Research  Range,  where  EXCEDE:  Spectral  was  conducted.) 

The  method  adopted  for  correcting  the  nominally  radially- 
symmetric  volume  emission  rates  for  these  off-perpendicular 
projections  is  described  later. 

Exposure  of  the  negative  is  within  its  dynamic  range, 
and  the  photographs  are  otherwise  of  good  quality,  at  all 
ejection  altitudes.  Characteristic  frames  are  reproduced  in 
montages  in  Ref's  1  and  13.  The  images  show  that  sensibly 
all  the  emission  in  the  field  of  view,  which  intercepts  the 
beam  axis  1-1/4  m  out  from  the  rocket  body,  comes  from  within 
-2  m  from  the  injection  flux  tube.  That  is,  no  general  dc 
corona  glow  surrounding  the  rocket  and  outside  the  beam  area 
proper,  such  as  might  be  ignited  by  energized  ionospheric 
electrons  returning  to  charge-neutralize  the  vehicle  (Ref  15), 
is  detectable.  Analysis  of  the  complementary  groundbased 
images  for  -125  km  rocket  altitude  showed  that  the  trans¬ 
verse-integrated  radiance  along  the  field-confined  beam  also 
lies  above  predictions  of  an  independent-particle  energy- 


dissipation  model  between  -40  m  (their  nominal  spatial 
resolution)  and  400  m  from  the  accelerator  (Ref  3).  On  the 
other  hand  the  end-point  range  of  the  beam  does  not  appear 
substantially  different  from  what  is  expected  (although  this 
issue  has  not  been  addressed  quantitatively;  refer  to  Ref  17), 
which  would  indicate  that  no  major  discrete  change  in  the 
energy  distribution  beam  of  the  electrons  can  be  occurring 
near  the  accelerator. 

When  the  rocket  is  above  -100  km  the  negatives  show  a 
broad  bright  area  extending  about  halfway  up  the  frame ,  from 
which  extends  an  apparently-narrowe r  region  whose  longitudinal 
increase  in  brightness  is  at  least  qualitatively  consistent 
with  the  increasing  sight  path.  In  retrospect  this  wider  area 
can  be  seen  to  be  associated  with  the  unexpectedly  steep  in¬ 
crease  in  glow  brightness  toward  the  rocket  within  the  first 
-5  m,  which  is  shown  in  the  isophote  plots  and  Fig  16  of  Ref  2. 
As  we  will  see,  the  broadening  turns  out  to  be  a  photographic¬ 
thresholding  and  perspective  effect,  since  the  halfwidths  of 
surface  brightness  and  volume  emission  rate  in  the  glow 
actually  decrease  toward  the  rocket  (closer  than  -10m  from 
the  accelerator).  The  lack  of  clear  dependence  on  ambient 
density  of  the  axial  radiances  in  this  close-in  region  over 
a  two  order  of  magnitude  range  in  density  is  one  of  the  more 
noteworthy  aspects  of  the  data;  the  emission  is  by  and  large 
consistent  with  some  near-constant  fixed  fraction  of  the 
beam's  kinetic  energy  being  absorbed  in  this  region. 

SELECTION  AND  INITIAL  REDUCTION  OF  THE  DATA 

We  selected  19  of  the  -270  useful  black-and-white  film 
images  for  analysis,  on  the  basis  of  the  following  criteria. 

The  set  includes  the  five  frames  previously  manually  reduced 
to  equi-radiance  contour  plots  and  reported  in  Ref  2;  these 
were  re-run  ab  initio  to  ensure  that  the  same  standard  micro- 


densitometry  procedures  were  applied  to  all  data. 

1)  Constant  ejected  current  (and  voltage)  from  a 
single  electron  accelerator  unit.  This  choice 
eliminates  dependence  of  the  photographic  radia¬ 
tion  yields  on  beam  energy,  and  furthermore 
maintains  a  fixed  viewing  geometry.  (EXCEDE's 
four  accelerators  were  spaced  50  cm  apart  in  a 
line  at  the  same  rocket  azimuth  as  the  cameras.) 

In  practice  unit  #4,  most  favorably  sited  farthest 
from  the  lens,  injected  7  to  8  amperes  over  most  of 
the  flight;  the  other  units  operated  sporadically. 
Some  information  on  the  volume  excitation  rates' 
dependence  on  beam  current,  which  although  appar¬ 
ently  not  needed  for  interpreting  EXCEDE:  Spectral's 
infrared  radiometry  data  would  be  useful  for  quali¬ 
fying  future  high-current  injection  concepts,  can 

be  derived  from  analysis  of  other  of  the  photo¬ 
graphic  frames. 

2)  Pairs  of  frames  at  the  same  altitude  (to  within 
~40%  in  ambient  air  density)  on  upleg  and 
downleg.  This  choice  serves  to  investigate 

the  effects  of  the  ejected  beam's  dwell  time 
on  individual  flux  tubes  (which  depends  on  the 
rocket's  velocity  relative  to  the  irradiated 
volume)  and  outgassing  from  the  vehicle  (which 
would  be  expected  to  decrease  with  time  after 
launch) . 

3)  Ejection-altitude  increment  less  than  a  scale 
height  of  N2»  -7-8  km.  Significant  changes 

in  the  energy  distribution  pattern  might  be  ex¬ 
pected  over  factor-e  changes  in  the  ambient  gas 
density. 


4)  Constant  film  exposure  duration,  78  millisec. 

This  limitation  removes  any  variability  due  to 
the  camera's  integration  time.  We  note  in  passing 
that  the  correction  for  reciprocity  failure  of 
the  film  (the  difference  in  effective  exposure 
among  the  8,  78,  and  800  msec  frames)  introduces 
only  negligible  error. 

Besides  not  considering  the  radiation  pattern  from  injection 
of  21-24  amperes  (3  accelerator  units  in  operation),  the  data 
set  selected  does  not  address  the  questions  of  variations 
between  adjacent  short-exposure  frames  or  long  and  short  ex¬ 
posures.  As  the  aforementioned  photometer  data  show  no  erratic 
10's  millisec-scale  variability  when  the  accelerator  current 
and  voltage  remain  constant  over  the  2-4  sec  ejection  pulses, 
reducing  adjoining  frames  would  determine  the  precision  of  our 
radiometric  photography.  The  fact  that  the  only  discernible 
change  in  the  photometer  signal  over  these  injection  periods 
can  be  ascribed  to  increasing  or  decreasing  air  density  leads 
to  the  expectation  that  the  spatial  radiance  distributions 
measured  in  long  and  short  film  exposures  would  be  similar. 

No  particular  effort  was  made  to  match  the  onboard  frames 
to  time  exposures  from  the  groundbased  cameras  (Ref  5). 

(These  provide  a  second  projection  of  the  glows  with  spatial 
resolution  much  poorer  close  to  the  rocket  but  far  better 
beyond  -40  m.)  This  is  because  the  limited  number  of 
high-quality  ground  photographs  available  would  restrict  the 
selection,  and  in  any  case  continuous  low  light  level  video 
images  exist  (Ref  17).  We  note  that  a  video  camera  was 
operated  onboard  EXCEDE:  Spectral,  which  due  to  what  appears 
to  be  its  restricted  dynamic  range  returned  very  little  useful 
information  (Ref  5). 


The  above  data-selection  conditions  were  met  for  alti¬ 
tudes  above  116  km,  where  on  upleg  accelerator  #4  achieved 
relatively  stable  operation.  An  additional  altitude -matched 
pair  results  from  inclusion  of  (one  of)  three  upleg  frames 
starting  at  106  km  during  which  the  same  unit  emitted  brief 
pulses.  These  images  provide  some  scaling  to  lower  injected 
currents,  but  have  poorer  signal/noise  because  of  their 
generally  lower  optical  densities.  (They  have  been  reduced 
but  not  yet  analyzed.)  Frames  meeting  conditions  1)  and  4) 
are  available  at  110,  98,  and  85  km  on  downleg;  these  are 
separated  by  almost  two,  rather  than  the  desired  less  than 
one  scale  height.  In  addition  we  included  two  8  millisec 
frames  at  79  and  74  km,  during  which  accelerators  #1  and  #3 
respectively  were  operating;  this  shorter  exposure  keeps  the 
film  densities  on  scale  at  the  higher  air  densities.  The 
anodes  of  these  two  accelerators  are  1-3/4  and  2-3/4  ro  from 
the  lens,  and  the  camera's  field  extends  to  within  5/8  m 
of  #1  (refer  to  Fig  12  of  Ref  1).  (These  also  are  not  yet 
evaluated,  as  they  have  a  different  photogramme  trie  scale 
than  the  other  frames,  the  payload’s  orientation  has  changed, 
and  the  end-point  or  "practical"  range  of  3  keV  electrons  at 
this  altitude  is  only  about  2m.)  Rocket  altitudes  at  the 
19  frames  analyzed  are  shown  in  Figure  4,  and  identifying 
data  for  them  are  listed  in  Table  1 . 

The  roll  of  1 6  mm  EK2475  film  had  been  processed  in  a 
reel  tank  using  agitation  known  to  result  in  spatially-uniform 
development.  Transmission  density  distributions  in  these 
frames  and  of  the  calibration  step  tablets  that  had  been  pre¬ 
viously  exposed  onto  the  film  (as  described  in  Ref  1)  were  read 
onto  seven-track  digital  magnetic  tape  by  our  scanning  micro¬ 
densitometer,  with  a  1 50ym  square  aperture  sampling  in  a 
rectangular  raster  on  a  48  x  50ym  grid.  This  oversampling 
results  in  approximately  43,000  run-length  coded  density  points 
per  film  frame.  The  digitization  increment  is  about  10%  in 
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Figure  4.  Altitudes  of  EXCEDE:  Spectral  78  msec-exposure 
photographic  frames  selected  for  evaluation. 

The  electron  accelerator  ejects  the  first  7-8 
ampere  beam  at  116  km  on  upleg.  The  diagram 
at  upper  right  is  a  side  view,  to  scale,  of  the 
camera  and  other  instrument  fields  and  (Ref  20) 
model  beam  interaction  volume. 


Table  1 .  Identification  of  EXCEDE:  Spectral  photographic 
frames  selected  for  evaluation 


TAPE  CT50.20  1  2  Feb  1982  Wedge  F 


1.  16163-1 

153.2  sec 

122.6 

km 

upleg 

2.  Step  tablet  2,  3  scans  step 

19,  2  scans  each  steps 

3.  16159-3 

136.0  sec 

115.7 

km 

upleg 

4.  16184+3 

246.5  sec 

111.8 

km 

downleg 

5.  16189-1 

267.9  sec 

97.7 

km 

downleg 

TAPE  CT5010  17  Feb  1982  Wedge  F 

1.  16192+8 

aborted 

2.  16192+8 

283.3  sec 

84.9 

km 

downleg 

(acce lerator 

off  during  ■ 

3.  16164+9 

160  sec 

124.6 

km 

upleg 

4.  16177-4 

214  sec 

124.9 

km 

downleg 

5.  16169-3 

179  sec 

127.8 

km 

upleg 

6.  16173-3 

197  se  c 

127.8 

km 

downleg 

7. 


1  . 
2. 

3. 
“  • 

4. 

5. 

6. 


Step  tablet  1,  2  scans  each  of  steps  19  through  6 
TAPE  CT5070  26  March  1982 Wedge  F 


118.9  km  upleg 
122.2  km  downleg 
119.7  km  downleg 


16160+9  143  sec 

16179-1  223  sec 

16180+9  230  sec 

empty  file 
16181-1  faulty,  aborted 
16181-1  232  sec  118.8  km  downleg 

Step  tablet  2,  2  scans  each  of  steps  18  through  1 

TAPE  CT5000  29  March  1982  Wedge  F 


1.  16185-1 

250  sec 

109.8  km 

2.  16194-2 

289.8 

sec 

78.9  km 

3.  16195+2 

295.0 

sec 

73.7  km 

4.  Step  tablet  2, 

2  scans  each  i 

TAPE 

CT5060 

1.  16156+4 

124.2 

sec 

109.3  km 

2.  16155+8 

120.7 

sec 

107.2  km 

-.  empty  file 

3.  16155-4 

118.4 

sec 

105.7  km 

4.  Step  tablet  3, 

2  scans  each  i 

5.  16154+4 

1  15.8 

sec 

104.0  km 

downleg 


#1  only) 
'8  msec,  ; 
#3  only) 


31  March  1982 


pulse ,  wedge  D 
tort  accelerator 
pulse ,  wedge  D 

tort  accelerator 
pulse ,  wedge  D 


wedge  A.  Very  noisy,  no  calibration 


Explanatory  Information  for  Table  1 . 

Table  1  lists  the  frames  digitized  from  the  black  and 
white  16  mm  EXCEDE:  Spectral  film  onto  five  separate  magnetic 
tapes.  These  data  were  later  transferred  to  one  labeled  multi¬ 
file  data  set  on  half  inch  7-track  magnetic  tape  at  the  AFGL 
Computer  Center.  The  frame  identification  scheme  used  here  is 
derived  from  the  film  footage  numbers  appearing  on  the  film  edge 
As  these  footage  numbers  each  span  two  frames,  each  number  is 
allocated  to  that  frame  adjacent  to  the  least  significant  digits 
e.g.,  frame  #16163  identifies  the  frame  adjacent  to  the  digits 
63.  Frame  #16163+5  is  five  frames  from  that  position  in  the 
direction  of  increasing  footage  numbers,  and  Frame  #16163-3 
is  three  frames  from  16163  in  the  direction  of  decreasing 
footage  numbers. 

All  the  photographic  step  tablets  used  for  calibration 
were  exposed  with  PhotoMetrics '  sensitometer  set  on  Stop  1024, 
through  a  Schott  BG-12  filter.  Tablet  1  is  a  1  sec  exposure, 

2  is  1/15  sec,  and  3  is  1/125  sec.  Wedges  F  and  D  are  contin¬ 
uous  density  reference  scales  incorporated  in  the  microdensi¬ 
tometer.  Their  density  quantization  steps  are  0.032  and  0.025 
units  respectively. 

The  data  on  the  labeled  multi-file  tape  are  stored  in 
their  original  form;  no  decoding  was  done  prior  to  recording 
on  this  tape.  Each  data  file,  i.e.,  each  digitized  image, 
is  stored  as  a  separate  file  with  a  unique  label.  All  files 
actually  containing  any  information  on  the  original  tapes 
were  loaded  onto  the  labeled  tape.  These  included  one  run 
aborted  during  digitization  (on  CT5010)  and  another  which 
was  believed  to  be  faulty  and  so  repeated  (on  CT5070) .  Any 
empty  files,  for  example  those  formed  by  consecutive  end-file 
marks  on  the  original  tapes,  were  removed  during  transfer  to 
the  labeled  tape;  these  empty  files  are  indicated  in  Table  1. 

The  labeled  tape  is  identified  by  the  number  OS0762  and 
the  system  of  labeling  individual  files  is  illustrated  by  the 
following  example, 

Label  EX264  is  derived  from: 

EXCEDE 

YE th  March  1982  --  digitization  date 

“4th  image  digitized  on  26  March  1982. 
Individual  data  images  may  be  accessed  and  copied  to  disk 
files  for  analysis  by  procedures  described'  in  the  Control  Data 
Corporation  NOS/BE  system  manual  under  the  general  heading  of 
labeled  multi-file  tapes. 


scene  radiance  over  most  of  the  image.  Frame  edges  were 
aligned  parallel  to  the  microdensitometer's  scan  direction 
to  within  one  scan  line  advance  by  reference  to  edges  of 
the  sprocket  holes.  These  holes,  besides  establishing  the 
angular  orientation  of  the  images,  also  were  used  to  deter¬ 
mine  the  frame  centers,  which  as  noted  align  with  the  optic 
axis  of  the  lens  in  the  long  dimension  of  the  frames  to  less 
than  one  scanning  slit  width. 

Software  was  written  for  reducing  this  film  density  file 
using  AFGL's  CDC-6600  digital  computer.  Programs  access  and 
identify  the  densities,  convert  them  to  scene  radiances  with 
the  view  angle -dependent  correction  for  lens  vignetting  (this 
had  been  done  manually  for  the  isophote  plots  reported  in 
Ref  2) ,  and  determine  axial  and  radial  distributions  of 
volume  emission  rate  using  the  unfold  method  described  below. 
Approximately  10  sec  of  central  processor  time  is  required 
to  decode  the  data  file  to  access  one  photographic  frame, 

2  sec  to  locate  its  center,  and  30  sec  to  compute  a  set  of 
volume  emission  rate  distributions  in  ten  cuts  through  the 
glow  transverse  to  its  long  dimension. 

Fig's  5  and  6  show  2-D  plots  of  raw  film  density  and  scene 
radiance  before  and  after  the  vignetting  correction,  for  a 
typical  frame.  Only  one-ninth  the  number  of  data  points  actu¬ 
ally  used  in  the  film  evaluation  are  retained  for  this  display. 
Note  (Fig  6b)  that  straight  lines  perpendicular  to  the  beam 
axis  and  the  view  direction  generally  map  onto  the  image  plane 
as  curved  lines;  in  addition  equally  spaced  sampling  points  in 
perpendicular  cuts  are  unequally  spaced  along  these  lines. 

As  the  viewing  directions  computed  for  individual  steps  in  the 
transverse  scans  do  not  in  general  correspond  exactly  to  posi¬ 
tions  of  digitized  data  points  in  the  film  plane  ,  we  devised  a 
simple  interpolation  scheme  to  determine  radiance  at  the  desired 
image  position  from  its  four  immediately-surrounding  stored  data 


(b) 

Figure  5.  Equi-density  contour  plot  (b)  and  print  (a) 

from  frame  at  127.8  km  on  upleg  (identification 
16169-3).  The  least  significant  digit  of  the 
optical  density  code  (integers  0-63)  is  printed. 
The  astigmatism  in  the  computer  plot  results 
only  from  the  character  spacing  of  the  printer 
and  in  no  way  affects  the  data  analysis. 
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(a) 


(b) 


Figure  6.  Uncorrected  (a)  and  corrected  (b)  scene  radi¬ 
ance  plots  from  frame  at  127.8  km  on  upleg 
(16169-3).  The  least  significant  digit  of 
the  radiance  code  is  printed.  The  ten  lines 
indicate  projections  perpendicular  to  the 
geomagnetic  field-aligned  beam  axis. 


DISTANCE  FROM  ACCELERATOR, 


points.  This  procedure  also  also  smooths  the  quantization 
noise  inherent  in  the  original  digitization. 

UNFOLDS 

The  cylindrical  symmetry  and  optical  thinness  of  the 
radiation  patterns  allows  their  volume  emission  rates  to  be 
unfolded  from  their  surface  radiance  distributions.  Column 
emission  rate  B(x)  at  points  intercepting  the  symmetry  axis  at 
transverse  distance  x  --  the  measured  quantity  --  is  related 
to  volume  emission  rate  V(r)  at  radial  distance  r  from  the 
symmetry  axis  by 

B(x)  =  2 fQr'  V(r)  r (r2  -  x2)‘1/2dr, 
where  V  goes  to  zero  at  r  =  r".  This  is  Abel's  integral,  and 
inverts  to 

V(r)  =  ir“Vrr'  (dB/dx)  (x2  -  r2)_1/2dx. 

We  extracted  linear  arrays  of  column  emission  rate  B(x)  at  ten 
ranges  D  along  the  beam  axis  from  the  electron  accelerator 
from  each  each  two-dimensional  image  brightness  array  (such 
as  Fig  6b),  increasing  x  in  equal  increments  in  object  space. 

We  first  attempted  to  fit  B  with  sets  of  third-  or  fourth- 
order  polynomials,  which  allows  Abel's  integral  to  be  solved 
analytically  (see,  for  example,  Ref  18).  The  line  data  sets 
were  organized  into  seven-point  segments  with  a  three-point 
overlap  at  each  end,  for  which  coefficients,  of  xm  were  calcu¬ 
lated  by  a  least-squares  method.  This  standard  procedure  was 
found  unable  to  smooth  the  occasional  large  (~20%)  fluctua¬ 
tions  in  radiance,  as  was  evidenced  by  persisting  and  unsys¬ 
tematic  spatial  structure  in  the  unfolds.  This  appears  to 
be  because  polynomials  tend  to  follow  large  excursions. 

Little  if  any  improvement  was  effected  by  varying  the  segment 
length  or  by  prior  smoothing  using  three-point  weighted  average 

As  the  transverse  radiance  distributions  look  to  be  roughly 
gaussian  in  shape  --  see,  for  example,  Fig  21  of  Ref  1  --  we 


substituted  fitting  by  single  gaussian  functions.  This  turned 
out  to  result  in  reasonable  fits  to  the  experimental  distri¬ 
butions,  as  evidenced  by  chi-squared  tests  (as  prescribed  in 
Ref  19  (see  also  Section  2);  the  chi-squared  or  residual  for 
unweighted  data  is  calculated  in  conjunction  with  each  least- 
squares  fit).  Figure  7  shows  a  computer  fit  to  one  of  the 
measured  line  radiances  in  Fig  6b.  Note  the  two  errant  data 
points  at  high  film  exposure,  which  we  believe  to  be  the  source 
of  the  unsuitability  of  the  more  general  polynomial  fitting 
method.  A  gaussian  offers  the  further  advantage  of  character- 
rizing  the  distributions  with  two  readily-recognizable  descrip¬ 
tors,  peak  surface  radiance  B0  (or  volume  emission  rate  VQ) 
and  width  parameter  a. 

Straightforward  unfold  of  a  gaussian  radiance  having 
full  width  at  half-maximum  2.354 a  gives  the  result  that  the 
volume  emission  rate  at  radius  r  is 

V  =  (2*a2)“1/2  B  =  (2w02) /2  Bq  exp-(r2/2a2) . 

The  assumed-constant  background  in  the  radiance  pattern 
(see  Fig  7)  is  removed  by  our  fitting  and  unfold  procedure. 

Note  that  the  brightness  and  volume  emission  patterns  have 
the  same  halfwidth  parameter,  which  as  we  have  defined  it 
makes  VQ  =  B0/CT/2ir  and  places  64%  of  the  total  emission 
within  distance  /2a  from  the  symmetry  axis,  where  the 
volume  emission  rate  is  V0/e. 

Gaussian-shaped  oblique  column  emission  projections  can 
not  of  course  result  from  strictly  gaussian  radial  volume 
emission  rate  distributions  whose  VQ  and  a  change  non- 
linearly  with  axial  distance  D  from  the  accelerator.  None¬ 
theless  insofar  as  the  fractional  changes  in  VQ  and  a  are 
not  large  over  individual  slant  paths  (column  emission  in 
these  paths  is  in  any  case  weighted  to  the  region  within  la 
of  the  glow's  axis),  the  radial  volume  emission  rates  can  be 
fit  to  gaussians  with  only  somewhat  larger  chi-squareds . 


PEAK  VOLUME  EMISSION  RATE  V0  ,  PHOTONS  M  SEC 


• 127.8km  Upleg,  69-3 
a,"  118.8km  Downleg,  81-1 

open  symbols:  uncorrected  values 
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Figure  8.  Examples  of  the  effect  on  V0 and  o  of  the  first- 
order  view  angle  correction  for  cases  where  the 
correction  is  small  (127.8  km  upleg)  and  large 
(118.8  km  downleg). 


«vy.  < 


We  made  a  first-approximation  correction  for  the  oblique  views 
using  the  following  procedure.  We  summed  three-dimensional 
volume  emission  rates  from  the  set  of  zeroth-order  VQ ,  a  deter¬ 
mined  for  D  =  1-1/2,  2,  3,  4,  5,  6,  8,  10,  15,  and  20  m  to 
derive  synthetic  transverse  radiance  distributions  viewing 
both  perpendicular  to  the  glow  axis  and  from  the  camera’s 
actual  position.  These  were  then  fit  to  gaussians  to  deter¬ 
mine  new  V0's  and  a's  for  the  two  sets  of  projections. 

The  improved  VG ' s  and  o' s  are  then  each  taken  as 

measured  zeroth-order  Vn  or  a 

X 

VQ  or  a  derived  from  synthetic  slant-path  data 
(VQ  or  a  derived  from  synthetic  perpendicular-projection  data) 
In  practice,  this  first  iteration  resulted  in  only  a  small 
change  in  VQ,  a  up  to  D  =  8  m,  and  in  some  of  the  frames 
a  large  (factor  -2)  decrease  in  a  and  increase  in  VQ  at 
D  =  15  and  20  m.  Cases  with  small  and  large  corrections  are 
illustrated  in  Figure  8.  A  second  iteration  may  be  indicated 
for  these  latter  data. 


RESULTS 

The  first-approximation  corrected  VQ ' s  and  a's  for  the 
four  upleg-downleg  pairs  are  plotted  in  Figures  9-12,  and 
plots  for  individual  downleg  frames  at  the  -lower  altitudes  of 
110  and  85  km  are  in  Figures  13  and  14.  Results  for  125  km 
downleg  (Fig  10)  are  qualitatively  different  from  the  others 
(a  is  larger  and  VQ  substantially  smaller  where  D  >  5m), 
which  suggests  that  this  data  record  should  be  re-examined 
for  unanticipated  errors.  The  altitude  dependence  of  VQ  and 
a  at  three  distances  along  the  beam  axis  from  the  accelerator 
is  shown  in  Figure  15. 

At  the  higher  altitudes  VQ  has  two  almost  straight-line 
segments  on  the  log-log  scale  and  is  largely  independent  of 
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127.8  km  Upleg,  69-3 
*,■  127.8  km  Downleg,  73-3 


(CORRECTION  SMALL) 
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Figure  9.  Photographic  volume  emission  rate  V0  on  the 
beam  axis  and  gaussian  width  parameter  o, 
for  altitude  127.8  km  on  upleg  and  downleg. 
The  sight  path  correction  is  small,  as  shown 
in  Fig  8. 
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Figure  10.  Photographic  volume  emission  rate  V 
on  the  beam  axis  and  gaussian  widthu 
parameter  a,  for  altitude  124.6  km 
upleg  and  124.9  km  downleg. 
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Figure  14.  Photographic  volume  emission  rate  V  on 
the  beam  axis  and  gaussian  width  para¬ 
meter  a,  for  altitude  84.9  km  on  downleg 


altitude,  while  o  at  constant  D  decreases  slowly  with 
altitude.  For  D  <  ~4  m  the  logarithmic  slope  of  VD  is 
about  -4,  and  beyond  ~5  m  it  falls  to  about  -f).3  (and  is 
less  well  defined).  The  distance  to  the  intercept  of  the 
two  manually-drawn  slopes  shows  an  increase  with  altitude, 
being  3,  3,  4,  5,  5  at  110,  118.  112,  125,  128  km  respec¬ 
tively.  No  sharp  change  in  slope  is  seen  at  85  km,  where 
the  practical  range  of  3  kV  electrons  is  only  12  m.  The 
large  difference  between  the  two  dependences  of  VG  on 
distance  to  the  accelerator  anode  suggests  that  different 
physical  processes  are  occurring  in  the  two  spatial  regimes. 

Even  in  the  low-slope  region  VQ  is  substantially  greater 
than  predicted  by  single -particle  energy  dissipation  models 
(Ref  3)  above  at  least  110  km.  In  any  case  the  theory  shows 
essentially  constant  radiance  viewing  perpendicular  to  the 
beam  direction  to  well  beyond  20  m. 

A  plausible  preliminary  interpretation  of  these  findings 
is  that  a  discharge  is  occurring  in  the  in-close  region, 
from  which  energetic  electrons  flow  down  the  field  lines  to 
produce  the  excess  excitation  beyond  about  3-5  m.  We  return 
to  this  interpretation  after  reviewing  further  features  of 
the  volume  emission  rate  geometry. 

The  gaussian  width  parameter  generally  starts  to  decrease 
with  distance  along  the  beam  at  D  =  6-10  m,  the  exception  being 
the  127.8  km  upleg  frame  (Fig  9).  While  the  o's  are  re¬ 
liable  when  the  view  angle  to  the  beam  axis  is  >45°  ,  when  D 
becomes  as  large  as  10  m  some  of  the  oblique-view  corrections 
are  quite  large  and  the  treatment  would  benefit  from  a  further 
iteration.  Nonetheless  the  turn-down  in  0(D)  is  present 
even  before  the  correction  is  applied  and  does  appear  to  be 
real.  It  may  be  caused  by  refocusing  of  the  initially  diverg¬ 
ing  beam  by  the  near-parallel  magnetic  field  as  noted,  expected 
to  occur  at  D  =  21  m.  This  would  mean  that  the  relative 


Figure  15.  Altitude  profiles  of  Vc  and  a  at  distances  2,  4  and  10  m  from  the 
accelerator.  The  dashed  lines  join  values  obtained  for  upleg  an 
downleg  at  the  same  altitude.  Except  at  125  km  the  upleg  Vc  and 
are  always  greater  than  on  downleg. 


velocity  spread  of  the  electrons  is  not  large,  that  is,  that 
collective  interaction  results  in  only  moderate  pitch  angle 
scattering  and  heating  of  the  beam.  On  the  other  hand  as 
the  primary  electrons  themselves  produce  only  a  fraction  of 
the  total  excitation,  their  crossing-over  would  not  be  ex¬ 
pected  to  have  so  substantial  an  effect  on  the  measured  spread 
parameter.  The  total  excitation  per  unit  path 
(/V  x  d  (area)  =  2 VC(D),  not  plotted)  is  decreasing  with 
increasing  D  beyond  6  m,  as  might  be  expected. 

o  as  well  as  V0(D  >~5  m)  can  be  seen  to  be  higher  on 
upleg  than  downleg,  the  exception  being  the  suspect  data  frame 
at  125  km.  This  effect  was  also  observed  by  the  onboard  photo¬ 
meter  (Ref  21).  The  volume  emission  rate  distribution  at 
85  km  is  qualitatively  different  from  that  at  higher  altitude, 
as  a  result  of  the  short  electron  range;  the  frame  at  97.7 
km  would  be  expected  to  show  a  transitional  behavior.  At  74  km 
lateral  Coulomb  scattering  is  expected  to  form  the  excitation 
into  a  near-spherical  volume  of  diameter  about  2  m,  which  with 
its  origin  from  accelerators  closer  to  the  camera  requires 
reprogramming  of  the  projection  geometry. 

COMPARISON  WITH  PRELIMINARY  MODEL  RESULTS 

Densities  and  temperatures  of  the  free  electrons  pro¬ 
duced  by  collective  interaction  of  the  ejected  charged  beam 
with  the  plasma  have  very  recently  been  computed  numerically 
(Ref  20).  Space -averaged  volume  rates  of  excitation  of  the 
major  electronic  and  vibrational  states  of  N£  were  calculated 
from  the  number  fluxes  and  energies  of  both  the  heated  secondary 
electrons  and  the  primaries.  The  model  assumes  that  a  beam- 
plasma  interaction  takes  place,  and  adopts  a  fixed  coupling 
coefficient  (0.1)  for  heating  of  the  plasma  electrons  (second¬ 
aries  and  ambient-ionosphere)  by  absorption  of  wave  energy 
from  the  resulting  rf  electric  fields.  To  date  it  averages 


this  energy  input  over  a  breakdown  oscillator  length  extending 
10  meters  along  the  geomagnetic  field  from  the  accelerator  and 
width  twice  the  Larmor  radius  of  3  keV  electrons  at  the  edge 
of  the  15°  half-angle  ejection  cone  (1.2  m)  ,  and  neglects 
transport  of  energy  out  of  this  volume. 

The  accuracy  of  the  preliminary  calculations  is  known  to 
be  limited  because  of  these  simplifications,  the  assumption 
that  the  electrons  have  a  Maxwellian  energy  distribution, 
and  omission  of  the  backscattered  and  hard  forward-scattered 
secondary  electron  component  (which  produces  about  as  much 
ionization  and  excitation  as  the  primaries  themselves  (Ref  10)). 
The  plasma  electrons'  velocities  would  be  expected  to  depart 
from  a  Boltzmann  distribution  near  the  higher  end  of  EXCEDE: 
Spectral's  altitude  range  because  the  time  between  electron- 
heavy  particle  (energy-loss)  collisions  becomes  comparable 
to  the  time  between  electron-electron  (energy-redistribution, 
and  thus  kinetic  tempe rature -establishing)  collisions.  (That 
is,  the  actual  population  reflects  directly  the  electrons' 
discrete  energy  losses  when  they  excite  ambient  atmospheric 
species.)  Nonetheless  the  model,  which  calculates  the  macro¬ 
scopic  results  rather  than  the  physical  origin  of  beam-plasma 
interactions,  is  a  source  of  ideas  for  interpreting  measured 
volume  emission  rate  distributions  such  as  in  Fig's  9-14  and 
emission  spectrums  in  terms  of  plasma-physics  processes.  The 
data  in  turn  serve  to  provide  input  parameters  --  such  as 
the  dimensions  of  the  interaction  region  --  for  refining 
the  calculations  of  electron  energy  distributions  and  thus 
excitation  rates  of  infrared  features. 

The  pertinent  output  of  the  calculation  is  time  depend¬ 
ent,  space -averaged  volume  excitation  rates  in  the  first 
10  msec  after  the  electron  beam  is  turned  on.  Although  its 
energy  deposition  is  of  course  spatially  continuous,  for 
purposes  of  comparing  theory  and  experiment  the  beam  may  be 
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Table  2.  Velocities  of  the  EXCEDE:  Spectral  rocket 
relative  to  the  geomagnetic  field. 


TIME 

TIHE-TOi 

ROCKET 

VELOCITY  (M/SEC) 

ALTITUDE 

PARALLEL 

PERPEND¬ 

TOTAL 

(SEC) 

(SEC) 

(KM) 

ICULAR 

136.0 

-51.5 

115.7 

399.1 

447.2 

599.4 

143.0 

-44.5 

118.9 

334.3 

433.4 

547.3 

153.2 

-34.3 

122.6 

240.3 

412.9 

477.7 

160.0 

-27.5 

124.6 

177.7 

399.2 

437.0 

179.0 

-8.5 

127.8 

3.6 

361.1 

361.1 

197.0 

9.5 

127.8 

-161.1 

325.7 

363.3 

214.0 

26.5 

124.9 

-316.9 

292.9 

431.5 

223.0 

35.5 

122.2 

-399.6 

275.5 

485.4 

230.0 

42.5 

119.7 

-463.8 

262.6 

533.0 

232.0 

44.5 

118.8 

-482.6 

258.1 

547.3 

246.5 

59.0 

111.8 

-616.8 

229.4 

658.1 

250.0 

62.5 

109.8 

-649.0 

223.0 

686.3 

267.9 

80.4 

97.7 

-816.9 

180.2 

836.6 

283.3 

95.8 

84.9 

-962.8 

128.8 

971.4 

289.8 

112.3 

78.9 

-1024.6 

99.8 

1029.5 

295.0 

107.5 

73.7 

-1070.2 

115.1 

1076.4 

i  APOGEE  IS  AT  TO  =  187.5  SEC  AFTER  LAUNCH 


visualized  as  entering  previously-und isturbed  air  volumes  in 
time  intervals  determined  by  the  interaction  volume's  dimen¬ 
sions  and  the  rocket’s  velocity.  For  orientation,  we  show  in 
Table  2  its  velocity  components  Vj^  and  vn  with  respect  to 
the  geomagnetic  field  direction  (refer  also  to  the  diagram 
in  Fig  4).  The  time  interval  in  which  the  beam  crosses  trans¬ 
verse  to  the  interaction  cylinder  is  -2.4/vj^  sec,  and  it 
passes  one  cylinder  length  down  or  up  the  field  in  ~10/vB  sec. 
(As  the  tabulated  velocities  show,  these  times  are  in  general 
comparable  to  the  10  msec  for  which  the  model  "freezes"  the 
beam  on  a  parcel  of  air.)  The  rocket  instruments'  fields  of 
view  as  well  as  the  beam  itself  dwell  longer  on  the  collective- 
interaction  volume  on  downleg  than  upleg:  the  fields  are  passing 
out  of  the  interaction  region  on  upleg  and  into  it  on  downleg 
(see  Fig  4),  a  point  to  which  we  shall  return  shortly. 

The  preliminary  model  results  for  EXCEDE's  electron 
ejection  conditions  show  short  (~3  msec) ,  intense  discharge 
pulses  near  apogee,  and  heating  of  the  secondary  electrons 
to  temperatures  near  6  eV  at  10  msec  at  the  lower  rocket 
altitudes.  This  temperature  decreases  by  about  5  eV  from 
its  narrow  peak  to  clamp  just  below  one  of  the  excitation 
potentials  of  N2  that  are  carried  in  the  code’s  tabulations 
of  energy-loss  cross-sections.  At  the  higher  altitudes  the 
rates  of  ionization  and  3914  S-band  emission  are  about 
a  factor  50  larger  than  those  due  to  individual  particle  im¬ 
pacts  during  this  pulse,  after  which  they  decrease  to  the  rates 
determined  by  single  particle  dynamics.  (As  noted,  the  back- 
flowing  and  hard  scattered  electron  component  is  not  considered 
in  the  model.)  The  tabulation  on  the  next  page  illustrates 
results  of  runs  to  date. 

The  interaction  model  in  effect  predicts  the  weak 
dependence  of  photographic  volume  emission  rate  on  ambient 
neutral-species  density  within  the  first  few  meters  from  the 


Iniection  altitude: 


Duration  of 
discharge  pulse 


Discharge  ionization/ 
beam  ionization,  max 

Plasma  electron  density, 
temperature  at  3  msec 

Average  3914  A  photon 
emission  rate ,  0-3  msec 


128  km 


~3  msec 


123  km 


~2  msec 


1  00  km 
No  avalanche 
deve  lops 


1 0x1 08/cm3 , 

21-1/2  eV 

1-1/2x1010/ 
cm3  sec 


5x108/cm3,  5x108/cm3, 


21  eV 


6.7  eV 


0.6x1010/  1-1/2x10^0/ 

cm3  sec  cm3  s pr 


cmJ  sec 


accelerator  (see  Fig  15).  (However  the  preliminary  absolute 
numbers  are  well  below  those  we  measured,  which  may  be  due 
to  numerical  errors  in  the  code's  plotting  routines.)  The 
existence  of  a  discharge  phase  with  duration  comparable  to 
the  dwell  of  the  beam  and  instrument  fields  on  parcels  of  air 
in  the  flux  tube  provides  a  potential  explanation  for  the 
differences  in  apparent  volume  emission  rate  between  upleg 
and  downleg  (as  noted  this  also  shows  in  the  onboard  3914  A- 
band  photometer  data).  Since  on  upleg  the  beam  passes  through 
undisturbed  air  volumes  at  a  higher  velocity,  and  the  fields 
of  view  move  more  rapidly  out  of  the  interaction  region,  the 
onboard  instruments  --  including  the  camera  --  see  on  the 
average  the  higher  air-excitation  rates  due  to  the  brief  dis¬ 
charge  pulse.  Note  that  the  increase  with  decreasing  rocket 
altitude  in  the  upleg/downleg  radiance  ratios  at  D  =  4  and 
10  m  (Fig  15),  and  also  in  the  photometer  data,  is  in  the 
same  direction  as  the  increase  in  ratios  of  the  vj/s  listed 
in  Table  2. 


Three  further  points  about  the  model's  predictions  merit 
notice.  First,  the  "zero-dimension"  calculations  do  not  of 
course  reproduce  the  steep  longitudinal  emission  gradients 
observed  in  the  first  10  meters.  (Space -re solved  calculations 
are  underway.  Ref  20.)  We  estimated  the  fraction  of  the  beam's 
energy  going  into  the  close-in  glow  by  extrapolating  the 
volume  emission  rate  distribution  at  128  km  to  D  <1-1/2  m 
from  the  accelereator .  The  result  indicates  that  at  least 
1/10  of  the  energy  --  the  figure  adopted  in  the  model  -- 
couples  to  the  plasma. 

Second,  even  when  no  discharge  is  ignited  the  secondary 
electrons  heated  by  plasma  instabilities  to  energies  above 
the  excitation  thresholds  of  electronic  and  vibrational  states 
of  ambient  and  newly-created  molecular  species  will  alter  the 
optical/ infrared  spectral  intensities.  Reference  is  made  to 
our  comments  in  Ref  1  on  this  issue,  in  particular  with  regard 
to  the  N2  triplet  sequence,  the  CO2  4.3ym  and  15um  bands 
(also  enhanced)  ,  and  the  unexpectedly  high  signals  from  water 
vapor  in  its  2.7um  and  6.3um  bands  and  perhaps  at  longer 
infrared  wavelengths.  We  point  out  that  the  measured  radiances 
in  these  and  other  infrared  emission  features  would  also  be 
different  at  the  same  altitudes  on  upleg  and  downleg  due  to 
the  effect  on  the  time -ave raged  signal  of  the  aforementioned 
temporal  variations  in  electron  temperature  during  the  beam- 
plasma  interaction.  (This  issue  has  not  been  directly  addressed 
in  the  spectrum  data  analyses  to  date.) 


LONGITUDINAL  EXTENT  OF  THE  EXCESS  EXCITATION 

Third,  these  heated  plasma  electrons  flow  out  of  the 
interaction  region,  in  both  directions  along  the  geomagnetic 
field  lines.  (Their  gyro  radii  are  much  smaller  than  their 
collision  mean  free  paths  down  to  at  least  85  km  altitude.) 
This  second  electron  source  enhances  the  Coulomb-backscatte r 


tail  behind  the  accelerator,  and  in  the  direction  in  which 
the  3  kV  beam  is  injected  it  both  increases  the  radiance  at 
wavelengths  where  photographic  film  and  video  cathodes  are  sensi¬ 
tive  and  --  more  important  for  EXCEDE  simulation  purposes  -- 
extends  the  distance  from  the  rocket  in  which  the  excess  exci¬ 
tation  referred  to  in  the  preceding  paragraph  takes  place.  The 
enhancement  in  photographic  signal  (over  that  predicted  by  an 
independent  particle  energy  dissipation-transport  model)  to 
a  few  hundred  m  at  128  km  injection  altitude  is  shown  in  Ref  3; 
note  that  the  cross-field  integrated  intensities  in  the  back¬ 
ward  direction  are  almost  as  large  as  those  in  the  ejection 
direction  (Fig  3-9)  ,  which  would  be  expected  if  high  fluxes 
of  ~21  eV  plasma  electrons  are  diffusing  outward  from 
the  interaction  volume. 

The  propagation  of  these  heated  electrons  can  of  course 
be  directly  computed  from  the  cross-sections  of  ambient  atmo¬ 
spheric  species  for  inelastic  collisions  and  elastic  scattering 
into  the  backward  hemisphere.  We  made  a  first-order  calculation 
assigning  energy  kT  to  each  electron  (that  is,  neglecting  the 
energy  spread  of  the  thermal  distribution)  and  assuming  that 
they  transport  in  a  one-dimensional  (field-confined)  random 
walk  due  to  >90°  collisions.  This  simplified  calculation  shows 
the  mean  free  paths  for  one  energy-loss  collision  to  be  1 20  m 
at  128  km  altitude,  60  m  at  123  km,  25  m  at  100  km,  and  5  m 
at  90  km.  The  bidirectional  flux  of  these  second-source 
electrons  reduces  to  that  of  the  unheated  secondaries  from 
the  beam  after  perhaps  three  such  interaction  lengths.  Thus 
their  effects  on  the  infrared  spectral  distribution  become 
small  beyond  D  >*75-350  m  above  100  km  beam  ejection  alti¬ 
tude,  and  at  substantially  shorter  distances  at  lower  altitudes. 

This  simplified  calculation  provides  estimates  of  the 
minimum  measurement-intercept  distance  for  remote  optical/ IR 
sensors  on  the  daughter  platform  in  mother-daughter  EXCEDE 


investigations.  The  electron  transport  calculation  can  be 
refined  and  even  specialized  to  excitation  of  specific  infra¬ 
red  emission  features  (for  example,  those  collisionally  ex¬ 
cited  by  very  slow  electrons)  by  applying  propagation  codes 
such  as  ARCTIC  (Ref  7).  This  effort  would  be  justified  when 
the  computations  of  the  plasma  electrons'  energy  spectrum, 
taking  into  account  the  dwell  time  of  the  beam  and  the  longi¬ 
tudinal  variation  in  energy  coupling,  are  finalized. 


A  NOTE  ON  OUTGASSING 

We  discuss  briefly  the  potential  role  of  outgassing  of  N2 
molecules  from  the  rocket  in  producing  the  measured  photo¬ 
graphic  volume  emission  rate  distributions.  That  the  radiation 
to  which  the  camera  is  responding  is  indeed  virtually  all  N2 
fluorescence  is  shown  by  the  glow's  spectrum,  and  secondarily 
by  the  measurement  of  3914  *  radiance  at  axis  intercept 
D  =  4-1/2  m.  To  result  in  the  observed  emission  pattern  the 
spatial  distribution  must  explain  (among  other  features)  the 
aforementioned  similar  decrease  of  V0  in  both  upleg  and  down- 
leg  ejection  geometries,  its  essentially  constant  near-in 
altitude  profile  above  100  km,  and  the  high  signals  from 
beyond  D  =  30  m  reported  in  Ref's  3  and  4.  Further,  the 
total  number  of  N2  molecules  calculated  as  outgassed  over  the 
trajectory  must  be  consistent  with  the  number  that  could 
reasonably  be  expected  to  be  carried  up  by  EXCEDE:  Spectral's 
0.6  m  radius,  6  m  long  Talos  Castor  vehicle. 

Desorption  of  much  more  easily  condensed  H2O  molecules 
from  EXCEDE  is  evidenced  by  the  strong  water  vapor  (and  some 
H2)  radiations  when  the  beam  was  on  (Ref  5) ,  and  perhaps  also 
by  collisional  quenching  of  other  emissions  by  H2O.  Its  y| 
and  v3  bands  at  2.73  and  2.66ym  (effective  radiative  life¬ 
time  10  msec)  w«.  ie  reported  detected  by  the  interferometric 
spectrometer  with  decreasing  intensity  between  at  least  105 


and  81  km  downleg  rocket  altitude  (p  133),  and  the  presence 
of  its  v2  band  at  6.3um  and  rotational  lines  above  17um 
(pp  183ff)  is  also  inferred.  (Preliminary  attribution  to  H2O 
of  ~2.75  -  3.0jjm  emission  seen  at  74  km  by  the  circular  vari¬ 
able  filter  spectrometer  (pp  29-30)  is  probably  in  error,  as 
the  unstructured  feature  lies  at  too-long  wavelengths.)  The 
interferometer’s  field  of  view  intercepts  the  beam  axis  at 
4  m  from  the  accelerator,  which  is  large  compared  to  the 
molecular  collision  mean  free  paths  in  this  altitude  range 
(see  below).  Thus  the  electron  beam  is  presumably  indirectly 
(by  excitation  transfer)  or  directly  (by  impact)  exciting  H2O 
left  in  the  rocket's  wake.  Altitude  profiles  of  intensities 
in  these  water  vapor  bands,  which  would  be  expected  to  differ 
on  upleg  and  downleg,  have  not  yet  been  derived  from  the 
spectrums . 

The  spatial  distribution  of  outgassed  molecules  calculated 
for  spherical-shaped  orbiting  satellites  (Ref's  22,  23;  other 
treatments  are  referenced  in  Ref  24)  can  be  adapted  for  esti¬ 
mating  the  dependence  of  outgassed  N2  concentrations  on  ambient 
air  density  and  distance  D"  from  the  cylindrical  sounding 
rocket.  The  model  predicts  that  the  density  of  desorbed  mole¬ 
cules  decreases  from  the  vehicle  surface  both  as  D'~( 1  or  2) 
and  exponentially  with  a  scale  length  about  0.4  times  the 
collision  mean  free  path.  This  length  is  550  cm  at  128  km 
altitude,  150  cm  at  120  km,  25  cm  at  110  km,  and  6  cm  at 
100  km  (and  1  cm  at  90  km  and  0.1  cm  at  80  km,  where  agreement 
of  volume  emission  rates  with  independent-electron  energy 
dissipation  calculations  has  not  been  tested  and  excitation  of 
water  vapor  is  still  observed).  These  numbers,  as  expected,  are 
comparable  to  the  initial-expansion  radii  of  luminous  chemical 
releases  (Ref  25).  As  the  rocket's  velocity  (Table  2)  is 
greater  than  the  molecules'  thermal  velocities  (after  colli¬ 
sion)  at  all  but  the  highest  ejection  altitudes,  the  close-in 
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spatial  distribution  is  different  in  its  wake  and  motion 
directions.  This  appears  to  contradict  the  observation  that 
the  volume  emission  rates  differ  between  upleg  and  downleg 
in  magnitude  rather  than  "shape"  (they  are  actually  lower 
looking  into  the  wake).  Further,  the  excess  high-altitude 
emission  at  D  *  5-20  m  when  the  electron  beam  is  directed 
either  ahead  of  or  behind  the  rocket,  and  also  at  !)  »  40- 
400  m,  falls  off  much  more  slowly  than  predicted  from  Ref  22's 
outgassing-transport  model. 

A  rough  estimate  of  the  outgassing  rate  required  can  be 
derived  by  comparing  an  N2  density  inferred  from  the  model 
with  that  which  would  fit  the  prediction  of  electron  energy 
deposition  theory.  At  123  km  and  D  (*D')  =  2  m  the  measured 
VQ  is  lOOx  that  predicted,  which  would  call  for  an  outgassed 
N2  density  1 0Ox  ambient,  or  3  x  10^/cm^.  Taking  either  Eq  12 
or  18  of  Ref  22  (for  spherical  spacecraft)  as  applicable  to 
EXCEDE,  for  which  we  adopt  an  effective  radius  of  1  m,  we 
find  that  the  required  outgassing  rate  is  about  1  mole/sec. 
This  means  that  to  fit  this  reasonably  typical  emission  data 
point  the  payload  and  recovery  system  must  be  releasing  as 
much  as  20  liters  of  N2  at  standard  temperature  and  pressure 
each  second. 

These  arguments,  while  by  no  means  the  last  word  on  the 
issue,  indicate  that  outgassing  alone  does  not  fully  explain 
the  spatial  and  altitude  profiles  of  N2  fluorescence  excited 
at  EXCEDE:  Spectral.  Note,  however,  that  outgassing  is 
thought  to  play  a  part  in  ignition  of  discharges  near  electron 
injection  rockets  (by  providing  a  more  dense  plasma,  Ref  26); 
thus  quantitative  calculations  of  the  desorbed  molecules' 
spatial  distribution  would  apply  in  refining  the  beam-plasma 
interaction  model. 


SUMMARY,  CONCLUSIONS,  RECOMMENDATIONS 

The  volume  emission  rate  distributions  of  N2+  and  N2 
fluorescence  unfolded  from  onboard  radiometric  photographs 
of  the  glows  produced  when  EXCEDE:  Spectral's  accelerator 
#4  injected  7-8  amperes  over  the  altitude  range  116  km  - 
128  km  apogee  -  85  km  show  the  following  principal  properties 

1)  Very  rapid  falloff  along  the  beam  axis 
within  the  first  ~5  m  from  the  rocket, 
with 

2)  little  dependence  on  ambient  density 
above  110  km  (and  perhaps  lower 
altitudes) ,  and 

3)  more  excitation  than  predicted  by 
individual-electron  energy  dissipation, 
theory  both  in  this  region  and  in  a 
region  of  slower  decrease  that  extends 
to  at  least  20  m; 

4)  a  fit  to  gaussians  transverse  to  the 
symmetry  axis,  with  width  parameter  that 
increases  up  to  -6  m  from  the  rocket 
and  then  decreases  in  all  but  one  or  two 
of  the  78  mill isec-exposure  frames  reduced; 
and 

5)  higher  maximum  values  and  width  para¬ 
meters  on  upleg  (that  is,  earlier  in 
the  flight)  than  on  downleg,  where  the 
electron  beam  and  fields  of  the  camera 
and  other  instruments  are  directed  into  or 
very  close  to  the  rocket's  wake. 

The  altitude  dependence  (Item  2)  and  existence  of  two 
emission  regimes  (Items  1,3)  are  consistent  with  a  discharge 
close-in  to  the  rocket,  the  excitation  in  which  is  being 
predicted  by  a  current  theoretical  study  program  (Ref  20). 


An  approximate  calculation  of  the  ranges  of  the  secondary 
electrons  that  flow  outward  from  this  region  in  both  direc¬ 
tions  along  the  geomagnetic  field,  using  the  mean  energies 
from  the  preliminary  model  calculations,  is  in  qualitative 
agreement  with  the  few  hundred  m  derived  from  an  analysis 
(Ref  3)  of  photographs  of  the  visible  glow  taken  from  a 
ground  station.  The  predicted  distance  along  the  injection 
flux  tube  beyond  which  IR/  optical  measurements  from  daughter 
platforms  would  include  only  negligible  amounts  of  radiation 
excited  by  plasma  electrons  --  that  is,  where  the  injection 
of  keV  electrons  more  precisely  simulates  excitation  of  the 
atmosphere  by  particles  from  nuclear  explosions  --  decreases 
with  decreasing  altitude  through  the  E  region.  However  the 
altitude  resolution  in  the  set  of  photographic  frames  so  far 
evaluated  is  insufficient  to  verify  estimates  of  the  altitude 
at  which  this  standoff  distance  approaches  the  long  dimension 
of  current  sounding  rocket  payloads. 

These  transport  calculations  can  be  readily  improved  and 
specialized  to  individual  impact-  or  indirectly-excited 
infrared  features  (for  example  those  from  CO2  molecules)  when 
the  energy  spectrums  of  the  heated  secondary  electrons  re¬ 
sulting  from  collective  interaction  of  the  primary  beam 
electrons  is  finalized.  Further,  the  measured  longitudinal 
and  lateral  distributions  of  excitation  provide  input  para¬ 
meters  for  model  calculations.  For  example  the  preliminary 
runs  assume  the  plasma  heating  regions  to  be  uniform  within 
in  a  10  m  long,  1.2  m  (Larmor-)  radius  cylinder,  while  the 
data  indicate  axial  limits  of  3  m  at  110  km  injection  altitude 
increasing  to  5  m  at  128  km,  and  gaussian  widths  decreasing 
with  increasing  altitude.  The  observed  extremely  strong 
dependence  of  emission  on  distance  from  the  accelerator's 
anode  indicates  a  need  for  space -resolved  calculations  of 
the  excitation  of  N2  and  N2+  by  plasma  electrons. 
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The  volume  emission  rate  patterns  do  not  appear  ascrib- 
able  to  outgassing  from  the  payload,  as  both  physically- 
unreasonable  spatial  distributions  and  total  number  of 
desorbed  N2  molecules  would  be  required.  Nonetheless  the 
question  of  outgassing  merits  further  consideration,  as  it 
affects  both  the  initial  molecule  densities  in  the  beam's 
path  and  (through  its  strong  H2O  component)  the  ionization- 
excitation  potentials  to  be  input  to  beam-plasma  interaction 
calculations.  We  have  provisionally  interpreted  the  higher 
radiances  on  upleg  (Item  5)  as  an  effect  of  the  ~3  msec 
(calculated)  duration  of  the  avalanche  phase  of  the  collective 
interaction  within  the  image  field.  This  idea,  however,  has 
not  been  quantified.  One  of  its  ramifications  is  that  the 
relative  spectral  distributions  at  ultraviolet,  visible,  and 
infrared  wavelengths  will  also  differ  along  EXCEDE's  trajec¬ 
tory,  a  point  that  has  not  been  investigated  from  the  data 
in  existing  emission  spectruras.  Indeed,  the  uncertainties 
inherent  in  the  model  calculations  indicate  a  pressing  need 
for  evaluation  of  the  upleg  and  downleg  altitude  profiles  of 
radiance  in  features  whose  excitation  mechanisms  are  quanti¬ 
tatively  understood  (for  example  the  N2  Second  Positive 
and  other  triplet  systems) ,  for  refining  the  energy  distri¬ 
butions  of  the  secondary  electron  component  that  is  also 
exciting  less  well  understood  infrared  features. 

Interpretation  of  the  decrease  in  gaussian  width  in  the 
outer,  low-slope  emission  region  as  due  to  refocusing  of  the 
initially-diverging  3  kV  beams  by  the  geomagnetic  field  does 
not  appear  consistent  with  the  conclusion  that  these  plasma 
electrons  contribute  much  more  of  the  fluorescence  excitation 
than  the  beam  electrons.  Beam-plasma  interaction  theory  should 
provide  further  information  about  the  velocity  spread  of  the 
primary  electrons  (which  would  be  expected  to  be  related  to 
the  energy  transferred  to  the  plasma) ,  and  about  the  physical 


interpretation  of  the  gaussian  width  parameter's  longitudinal 
variation.  In  this  regard,  the  end-points  of  detectable  glow 
radiance  in  ground  photographs  (which  are  only  partially 
analyzed)  so  far  do  not  indicate  severe  energy  degradation 
of  EXCEDE:  Spectral's  beams. 

The  major  conclusion  of  this  evaluation  of  EXCEDE: 

Spectral  closeup  images  is  that  excess  fluxes  of  secondary 
electrons  accompany  the  primary  field-confined  beam  to  at 
least  20  m  at  injection  altitudes  above  ~100  km.  As  the 
expected  energy  distributions  of  these  electrons  (near- 
Maxwellian,  with  peaks  comparable  to  the  excitation  potential 
of  atmospheric  species,  Ref  20)  differ  from  that  of  the 
secondaries  produced  in  impact  ionization  by  fast  charged 
particles  (power-laws.  Ref  10),  the  infrared  emission  spec- 
trums  they  excite  in  the  atmosphere  will  differ  out  to 
10's  -  100's  m  from  the  rocket  at  altitudes  where  conventional 
thermionic  cathodes  and  cryocooled  spectroradiometers  operate 
without  severe  interference  from  the  atmosphere.  Comparison 
of  the  measured  fluorescent-emission  patterns  with  early  model 
predictions  is  still  preliminary,  and  several  of  the  geometric 
features  have  not  been  interpreted. 
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SECTION  2 


CORRELATION  BETWEEN  2. 8-3.1  vim  CHEMILUMINESCENCE 
AND  PARTICLE  ENERGY  INPUT  TO  THE  ATMOSPHERE 

BACKGROUND 

The  high  spatial-  and  temporal-resolution  aircraft  meas¬ 
urements  of  short  wavelength  infrared  emission  from  aurora 
(Ref's  1,27)  represent  the  primary  source  of  direct  informa¬ 
tion  about  this  radiation's  correlation  with  deposition  of 
ionizing  energy  in  the  real  atmosphere,  and  thus  about  struc¬ 
ture  of  the  SWIR  sky  background  resulting  from  nuclear  explos¬ 
ions.  In  this  section  we  quantify  the  relationship  between 
near-zenith  radiances  in  the  2.832-3.125  (FWHM)  wavelength 
band  (the  infrared  output,  from  vibrational  chemilumin¬ 
escence  of  NO)  and  a  fluorescence  bajid  of  N2+  (which  measures 
the  energy  input  rate  into  the  air  column)  using  the  procedures 
we  developed  as  part  of  the  Test  Plan  for  the  aircraft  data 
missions  (Ref  28). 

The  measurements  were  made  with  a  cryocooled  radio¬ 
meter  accurately  coaligned  (as  described  in  Ref  29)  with  a 
conventional  photoelectric  filter  photometer  having  the 
same  6  millirad-square  field  of  view.  (The  atmosphere  above 
the  NKC-135A  jet  aircraft's  operating  altitude  is  essentially 
transparent  and  presents  only  small  thermal  background  in 
the  selected  sensitivity  band.)  This  field  gives  a  7/10  km 
instantaneous  footrint  at  auroral-emission  altitudes,  compar¬ 
able  to  that  under  consideration  for  future  spaceborne  infrared 
surveillance  systems.  Signal  integration  time  for  S/N  -  10 
was  less  than  1  sec,  during  which  the  aircraft  moves  2/10  km. 
The  instruments  were  fixed-mounted  in  AFGL/DNA's  Optical 
Flying  Laboratory  (serial  55-3120)  to  point  to  the  magnetic 
zenith  when  it  was  flying  south  along  a  geomagnetic  meridian 


at  auroral  latitudes  in  the  Western  Hemisphere,  so  as  to  view 
particle  beams  restricted  to  a  narrow  flux  tube  and  thus  with 
a  nominally  unique  and  calculatable  altitude  profile  of  energy 
deposition.  Data  taken  when  the  aircraft  was  returning  north¬ 
ward  in  series  of  passes  underneath  the  auroral  oval,  and 
from  some  west-east  flight  trajectories ,  were  also  analyzed. 
Support  instrumentation  onboard  included  a  photographic  all¬ 
sky  camera  (images  spanning  the  useful  radiometer  data  seg¬ 
ments  are  shown  in  Ref  1),  all-sky  and  narrow-field  video 
camera  with  better  temporal  resolution  (so  far  not  applied 
in  the  data  evaluation) ,  hydroxyl  airglow  radiometer  (for 
measuring  the  comparably  intense  but  very  slowly  varying 
fundamental  Me ine 1-bands  background  in  the  2.8  -  3.1pm 
signal),  and  multi-channel  filter  photometer  (which  determines 
the  peak  energy  deposition  altitude  of  the  auroral  electrons 
(Ref  30) ,  and  in  addition  served  in  calibrating  the  primary 
air-fluorescence  photometer). 

Further  particulars  of  the  aircraft  measurements,  back¬ 
ground  information  on  the  aerochemistry  of  N0+ ,  and  program 
results  to  date  are  given  in  Ref's  1  and  27  and  the  references 
therein.  To  place  the  current  data  evaluation  in  context  we 
briefly  review  here  the  previous  findings  from  the  flights 
of  the  narrow-field  instruments  in  1979-1980  and  similar 
earlier  missions  with  a  160  mrad-field,  30  sec-response  pair 
whose  axis  pointed  to  the  aircraft's  zenith  (also  supported 
by  the  instrumentation  described  above).  The  yield  of  SWIR 
photons  in  the  radiometer's  passband  from  the  atmosphere 
between  150  and  100  km  is  linear  with  energy-input  rate  in 
the  range  reached  in  aurora,  ~105  -  1()6  io.n  pairs/cm^-sec ; 
however  this  energy  conversion  efficiency  varies  over  at 
least  a  factor  two  within  individual  flight  missions  and  a 
factor  3  among  missions,  for  reasons  that  have  not  been  ex¬ 
plained.  The  mean  yield,  and  even  the  minimum  yield,  is 


higher  than  predicted  from  laboratory  and  theoretical  inves¬ 
tigations  of  NO  vibrational  band  excitation/de-excitation 
(on  which  more  shortly) ,  but  consistent  with  the  results  of 
sidelooking  rocket  radiometry  of  auroral  forms  (Ref  30).  No 
statistically-signif icant  correlation  has  been  identified 
between  this  energy  conversion  efficiency  or  the  delay  in 
SWIR  emission  following  particle  input  and  geophysical  or 
instrument  factors  such  as  solar  illumination  of  the  auroral 
ionosphere,  rate  of  energy  input,  Kp,  measured  intensity  of 
the  hydroxyl  airglow  (which  might  indicate  calibration  in¬ 
stability)  ,  aircraft  flight  direction,  or  position  of  the 
flight  in  the  series.  On  the  other  hand  the  yield  of  NO 
photons  averaged  over  '-I -min  data  segments  appears  to  be 
increasing  with  altitude  of  the  maximum  in  the  auroral  electron 
energy  deposition  profile,  in  the  110  -  130  km  altitude  range 
containing  most  of  the  peaks'  occurrences  (Figure  16). 

The  aircraft  radiometry  data  have  resolution  far  superior 
to  that  of  the  SWIR  output-energy  input  measurements  from 
rockets  (or  contemplated  research  satellites,  which  have  less 
favorable  viewing  projections).  Quantification  of  the  temporal 
spatial  correlations  would  permit  the  visible -photographic  data 
on  past  high  altitude  nuclear  detonations  to  be  used  for  charac 
terizing  sky  clutter  in  the  important  2. 8-3.1  urn  wavelength  band 

DATA 

A  total  of  15  segments  from  10  flights  has  been  selected 
for  analysis  on  the  basis  of  their  large  infrared  signal  enhan¬ 
cements  coincident  with  visible  aurora.  The  ~70  min  of  data 
include  a  mix  of  instrument  pointing  directions  (5  segments 
are  from  the  intended  magnetic  zenith)  and  solar  depression 
angles.  Another  6  segments  (~15  min),  in  which  the  column 
emission  rates  in  the  N2+  (0,0)  First  Negative  fluorescence 
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Figure  16. 


AVERAGE  NO  OVERTONE  YIELD,  % 


Altitude  dependence  of  NO  overtone  emission 
yield  in  the  15  aicraft  data  segments.  The 
error  bars  on  yields  are  statistical,  and 
the  altitude  uncertainty  (+10,  -5  km)  is 
omitted. 
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band  ranged  from  20  up  to  50  kilorayle ighs ,  may  also  merit 
evaluation.  (Column  energy  input  is  proportional  to  intensi¬ 
ties  in  this  sequence  within  10%  over  the  measured  range  of 
peak  deposition  altitudes  (Ref  30;  see  also  Section  1).) 

The  signal/noise  and  signal/background  ratios  in  these  data 
segments  are  indicated  by  Fig  22 's  reproduction  of  the  typic¬ 
ally  1-sec  pre-integrated  and  sampled  traces;  we  applied  no 
further  filtering.  Flight  conditions  for  these  data  are 
listed  in  Table  3  (from  Ref  1). 

Chemi luminous  yields  in  the  nitric  oxide  overtone  system 
are  plotted  against  altitude  of  the  auroral  particles'  maximum 
energy  deposition  in  Fig  16.  Effective  height  of  the  columns 
is  about  30  km,  10  km  of  which  is  below  the  profile  peak.  The 
error  bars  indicate  90%  confidence  limits  on  the  slope  of  the 
best-fit  regression  line  between  simultaneously-measured  column 
energy  input  and  infrared  output  (as  in  Fig  4  of  Ref  27;  an 
averaging  method  was  used  as  a  check  in  cases  that  output 
lagged  input  by  more  than  ~3  sec)  ,  and  do  not  include  possible 
experiment  sources  of  error.  The  estimated  vertical  error 
bar  (two  standard  deviations,  not  shown)  is  +10,  -5  km. 

These  yields  assume  that  a  fixed  fraction  of  a  model 
spectral  distribution  in  the  nitric  oxide  overtone  system  -- 
42%,  principally  from  bands  originating  from  the  6th  to  12th 
vibrational  states  --  remains  in  the  wavelength  region  to 
which  the  radiometer  responds.  The  actual  spectrum  of  emission 
from  the  atmosphere  has  not  been  measured  with  sufficient 
resolution  to  identify  individual  &v  =  2  transitions.  It  is 
expected  to  be  weighted  to  these  measured  (longer-wave length) 
bands  at  lower  energy-deposition  altitudes  because  the  vibra¬ 
tional  cascade  is  partially  quenched  by  collisions  (Ref  31); 
however,  as  Fig  16  shows,  the  yields  in  the  instrument's  sensi¬ 
tivity  band  indicate  a  decrease  with  decreasing  altitude.  (We 
return  to  this  issue  in  the  following  subsection.)  The  obser- 
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vation  that  the  energy  conversion  efficiencies  that  were 
measured  with  narrow-field  instruments  are  generally  higher 
(by  almost  a  factor  two)  than  those  previously  reported  from 
the  wider-field  instruments,  like  the  spread  within  each  data 
set,  is  so  far  not  explained. 

The  output  of  SWIR  photons  is  found  to  lag  energy  input 
in  all  but  one  of  the  segments,  in  which  no  delay  is  detect¬ 
able;  in  no  case  is  the  SWIR  signal  seen  to  lead.  We  esti¬ 
mated  these  lags  (in  Ref's  1  and  27)  from  the  offsets  between 
discrete  radiance  peaks  and  from  differences  in  time  for  the 
two  signals  to  return  to  baseline  when  excitation  cuts  off, 
finding  that  they  vary  from  0  to  6  sec  with  the  distribution 
to  be  shown  later  in  Fig  19  (and  as  noted,  without  any  signi¬ 
ficant  correlation  with  geophysical  or  instrument  factors). 
Delays  interpreted  as  the  inverse  of  NO*  quenching  rates  were 
also  computed  (in  Ref  27,  by  another  group)  using  a  method 
that  involved  dividing  complex  Fourier  amplitude  transforms 
of  the  input  and  output  data  streams  to  derive  a  transfer 
function  of  the  atmosphere.  This  array  was  filtered  after 
retransforming  into  the  time  domain  applying  an  unspecified 
smoothing,  and  then  was  again  transformed  back  into  frequency 
space  and  fit  to  a  single-or  double-pole  filter  representation 
This  procedure  results  in  quenching  coefficients  all  but  two 
of  which  lie  within  1.1  ±  0.1 /sec.  (The  result  for  one 
segment,  for  which  we  had  estimated  0  to  1  sec  delay,  was 
2.0/sec.)  This  clustering  may  result  from  the  use  of  a 
smoothing  window  of  effective  width  0.9  sec.  Before  deter¬ 
mining  delay  times  by  applying  Ref  28 's  cross-correlation 
and  curve -fitting  methods,  we  discuss  briefly  the  existing 
information  about  excitation  and  quenching  of  NO* . 

EXPECTED  SPATIAL- TEMPORAL  CORRELATION 

The  current  model  for  production  of  vibrationally- 
excited  NO  molecules  following  ionization  and  dissociation  of 


N2  In  the  upper  atmosphere  is  illustrated  schematically  in 
Figure  17.  Oxygen  atoms  can  be  seen  to  play  an  important 
part  in  the  spectral  yield  and  grow-in/decay  time  of  the 
radiation,  which  is  still  unsatisfactorily  quantified. 
Specifically , 

--  NO(v)  is  deactivated  by  0,  with  partial 
rate  coefficients  for  NO(v')  not  yet 
reliably  determined; 

--  n2d,  the  principal  precursor  of  NO*  ,  is 
deactivated  by  0  with  a  rate  coefficient 
whose  measured  values  (Ref's  32,33)  differ 
by  a  factor  5  at  120  km;  and 

--  N^D  is  produced  when  N2+  reacts  with  0, 
but  not  when  charge  exchanges  with  O2. 

In  addition  the  probabilities  that  N  atoms  in  the  state  result 
from  the  key  dissociation  reactions  are  subject  to  error.  A 
recent  calculation  (Ref  31)  indicates  that  the  uncertainties 
in  the  rates  of  the  first  two  processes  above  and  in  the  atmo¬ 
sphere's  [01  profile  result  in  a  factor-2  uncertainty  in  the 
predicted  steady-state  volume  emission  rates  of  NO  overtone 
at  altitudes  where  auroral  particles  (and  fission  debris  of 
nuclear  explosions)  deposit  energy.  However  the  highest  derived 
column  emission  rates  are  still  almost  a  factor  two  below  the 
mean  of  all  the  aircraft  (and  aforementioned  rocket)  measure¬ 
ments.  Runs  of  DNA's  ARCTIC  aerochemistry  code  (Ref  34)  also 
resulted  in  similarly  low  overtone  yields  under  steady-state 
particle  irradiation. 

Characteristic  times  in  which  the  zenith  auroral  radi¬ 
ance  approaches  asymptotically  these  equilibrium  values  have 
also  been  calculated  by  ARCTIC  (Ref  34),  using  model  energy 
deposition  profiles.  They  are  1-1/2  and  1  sec  when  the 
maximum  of  the  profile  is  at  100  km  and  the  input  is  stepwise 
increased  and  decreased  respectively,  and  7  sec  when  this 


nitric  oxide  vibrational  radiations  by  energetic  particles. 
ie  minor  reactions  have  been  omitted.  The  principal  reactio 
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maximum  is  at  120  km  and  the  input  is  increased.  AFGL/DNA's 
"BRIM"  model  of  infrared  limb  brightness  structure  (Ref  35) 
makes  the  simplifying  assumption  that  the  final  process  of  the 
principal  reaction  paths,  N^D  +  O2  N0(v)  +  0  (whose  rate 
coefficient  is  5.2  x  10~12  cm^/sec)  is  the  rate-limiting  step 
in  NO  chemiluminescence.  This  directly  implies  that  buildup 
of  N0+  ions,  whose  characteristic  time  for  reacting  with 
electrons  is  several  sec,  is  not  a  major  factor  in  production 
of  N^D  atoms  compared  to  direct  impact  dissociation  of  N2 
molecules.  With  this  assumption  the  time  dependence  is 
exponential  with  time  constant  tq  =  (5.2  x  10“^ 
which  is  1/10  sec  at  100  km,  4  sec  at  120  km,  and  70  sec  at 
150  km  deposition  altitude.  The  infrared  volume  emission 
rate  at  time  t  would  then  be 

(steady-state  energy  conversion  efficiency)  X 

To*1  I(t')  exP  ’  Kc  -  t')/xD]dt', 

where  l(t')  is  the  energy  input  rate  at  previous  time  t'.  The 
vertical-column  emission  rate  would  be  this  volume  emission 
rate  weighted  over  the  auroral  deposition  altitude  profile. 

The  aircraft  measurements  --  and  for  that  matter  any 
passive  optical  remote  sensing  of  aurora,  except  perhaps  of 
narrow  isolated  arc  forms  --  of  course  refer  to  a  range  of 
emission  altitudes.  As  noted,  this  range  spans  about  30  km 
in  the  auroras  in  the  data  set;  it  decreases  with  decreasing 
altitude  of  the  profile  peak  (see,  for  example.  Ref  36)  and 
thus  generally  with  auroral  intensity,  which  statistically 
correlates  with  hardness  of  the  particle  energy  spectrum 
(Ref  37).  Thus  some  imprecision  is  expected  in  representa¬ 
tions  of  the  altitude -summed  time  delays  by  single  parameters 
such  as  we  derive  in  the  next  section. 

The  lag  attenuates  high-frequency  components  of  the 
temporal  frequency  spectrum  of  SWIR-band  emission  excited  by 
nuclear  particles,  with  a  transfer  function  that  can  be 


straightforwardly  calculated.  A  similar  effect  on  spatial 
frequencies  of  the  background  radiance  results  from  differ¬ 
ential  transport  of  precursor  species  by  wind  shear  and 
turbulent  diffusion  that  occurs  naturally  or  is  induced  by 
the  detonation.  The  horizontal  wind  at  120  km  has  a  shear 
that  is  highly  variable  in  magnitude  and  direction  (Ref  38) , 
which  averages  about  100  m/sec  per  10  km  vertically;  thus 
irradiated  air  parcels  at  different  altitudes  in  a  few  sec 
move  distances  relative  to  one  another  that  are  comparable 
to  the  projected  minimum  surveillance-sensor  footprint. 

DATA  REDUCTION  --  APPROACH 

To  make  quantitative  estimates  of  the  delays  we  apply 
two  methods  discussed  in  the  Test  Plan  (Ref  28) ,  established 
statistical  cross-correlation  and  a  procedure  developed  for 
determining  predosing  parameters  by  fitting  time-integrated 
power  input  to  observed  radiation  output  using  a  simplified 
model  of  the  latter's  grow-in  and  decay.  As  noted,  the  high 
signal/noise  data  segments  selected  are  identified  in  Table  3, 
whose  column  headings  are  described  in  Ref  1 ;  montages  of 
all-sky  views  of  the  visible  aurora  are  in  Ref’s  1  and  27. 

Cross-correlation  analysis  is  particularly  useful 
where  foreknowledge  of  the  physical  relationship  between 


processes  is  limited.  The  cross-correlation  function  is 
defined  (neglecting  for  the  moment  normalization)  by 
C(t )  =  I(t)  J(t  +  t ) d t , 

where  I(t)  is  the  input  (3914  A  photometer)  signal  at  time 
t,  and  J(t  +  t)  is  the  output  (2.8  -  3.1ym  radiometer) 
signal  at  later  time  t  +  t .  A  narrow  maximum  in  C  at 
t  *  x i ,  for  example,  would  mean  that  the  SWIR  emission 
is  shifted  by  tj  seconds;  in  practice,  we  expect  it  to  be 
delayed  over  a  range  of  times  comparable  to  that  required 
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for  the  processes  in  Fig  17  to  approach  completion.  Were 
J(t)  to  remain  proportional  to  I(t)  (no  delay),  C(t)  would 
be  directly  proportional  to  the  input's  autocorrelation 
function 

A(x)  =  f_J*  I(t)  I(t  +  x ) d t , 

which  has  a  maximum  at  t  =  0.  As  derived  in  Ref  28,  in 
the  more  general  case 

C(t)  =  (photon  yield)  X  f  T°°  R(t  +  t")  , 

where  R  is  a  predosing  function  to  be  discussed  shortly  and  t" 
is  the  time  interval  between  input  and  output. 

To  compute  C  and  A  we  manually  digitized  14  of  the  pairs 
of  data  traces  at  1  sec  intervals,  which  results  in  about  180 
points  per  segment;  segment  12  was  digitized  each  0.4  sec 
since  its  sharp  SW1R  enhancement  lasted  only  about  8  sec. 

Each  radiance  trace  was  then  time -averaged  to  zero,  and  a 
normalized  C(t)  calculated  as 

(P-1-0 

C(x)  =  (P  -  t)-1  [I  <t)  -  I]  [J  (t  +  x)  -  J]  At, 

t=0  _  _ 

Where  P  is  the  total  number  of  data  points  and  I  and  J  are 

the  average  values  in  the  interval  (the  factor  (P  -  x)~1  com¬ 
pensates  for  the  decreasing  number  of  terms  in  the  sum  as  lag 
t  increases).  Subtraction  of  I  and  J  formats  the  above 
equation  to  match  the  well  known  covariance  of  statistical 
analysis;  historical  usage  in  both  statistics  and  turbulence 
theory  applies  the  term  cross-correlation  to  the  normalized 
ratio  C(t)/C(0),  to  which  results  presented  in  the  follow¬ 
ing  subsection  refer.  We  computed  C(t)  for  x  =  0  to  10  sec  in 
1  sec  steps,  finding  that  almost  all  of  the  segments  displayed 
broadly- peaked  cross-correlation  functions  decreasing  back  to 
C(0)  in  less  than  10  sec. 

The  second  procedure  fits  parameters  in  a  model  that 
calculates  the  contribution  of  previous  energy  input  to  SWIR 


output.  It  is  summarized  in  the  equation 

qN0  =  F  R(t  " 

where  q^Q  and  qQ  are  volume  production  rates  of  NO*  and 
initial  ionizations,  F  is  the  fraction  of  these  ionizations 
that  result  in  NO  vibrational-band  photons,  and  R  is  an 
impulse -response  function  that  expresses  the  rate  at  which 
earlier  energy  input  leads  to  emission  from  NO.  (Transport 
of  precursors  is  neglected,  and  dissociation  of  N2  molecules 
is  lumped  with  impact  ionization,  to  which  it  is  in  any  case 
proportional.)  With  the  assumption  that  R  does  not  change 
substantially  over  the  ~30-km  range  of  energy  deposition 
altitudes  in  typical  aurora, 

qo(t')  =  K'I(t'),  qNO(t)  =  K'M(t), 

where  K"  and  K"  are  constants  (K""  is  presumed  to  depend 
on  the  column's  mean  altitude).  The  predicted  SWIR  column 
emission  is  then 

Jp(t)  =  Uj  KO  R(t  - 

where  K  lumps  the  proportionality  constants  K' ,  K" ,  and  F. 

This  equation  relating  the  data  streams  from  a  single 
coaligned  radiometer -photometer  pair  applies  to  aurora  that 
is  coherently  pulsating  while  remaining  spatially  near-uniform 
over  the  region  scanned  by  the  moving  aircraft  instruments. 

Large  rates  of  change  of  I  provide  the  highest  resolution  in 
unfolding  R,  and  constant  input  of  course  provides  no  informa¬ 
tion  about  predosing.  We  select  for  R  a  physically  reasonable 
analytical  form  with  two  free  parameters  that  are  fit  to  the 
data.  The  aerochemistry  model  (Fig  17)  indicates  that  follow¬ 
ing  impulsive  energy  input  R  increases  as  the  concentrations 
of  N^D  and  N0+  (and  perhaps  other  precursor  species)  build 
up,  and  then  decreases  as  these  species  become  depleted.  This 
suggests  a  function  of  the  form 

R  =  G(t  -  t')n  exp[-(t  -  t')/x0], 
where  tc  is  the  time  constant  of  the  decreasing  contribution 
of  earlier  dosing,  n  expresses  the  rate  of  grow-in  of  excitation 
of  NO* ,  and  G  normalizes  /®Q  Rd(t  -  t')  to  1.  (As  Ref  28 


shows,  G  is  a  gamma  function  and  reduces  to  1/tq  for  n  =  0.) 
The  equation  that  predicts  SWIR  radiance  then  becomes 
Jp(t)  =  KG^C  I(t')  (t  -  t')n  exp[-(t  -  t')/t0]. 

K,  which  refers  to  the  number  of  SWIR  photons  emitted  per 
3914  A  photon,  can  be  determined  directly  from  data  taken 
when  the  fractional  change  in  input  rate  is  small  over  the 
expected  predosing  time,  or  by  the  crossplot-averaging  methods 
of  Ref's  1  and  27.  Note  that  two  parameters  allow  independent 
rates  of  grow-in  and  die-out  of  NO  vibrational  emission;  as 
mentioned  BRIM  (Ref  35)  assumes  instantaneous  grow-in, 
with  n  =  0. 

To  evaluate  n  and  t0,  we  summed  the  above  integral  over 
a  t  -  20  sec  interval  (the  aircraft  platform  moves  4.4  km  or 
~5  instrument  fields  of  view  in  this  time).  The  process  was 
repeated  (by  computer)  for  n  and  td  =  0  to  10  in  one-unit 
steps,  with  smaller  steps  applied  near  the  best-fit  value. 
Best-fit  is  at  the  minimum  value  of  the  sum 

X2  =  M-1  I  fj(tt)  -  Jp(ti)]2/ai2, 

where  a*  is  the  uncertainty  in  the  measurement  of  J[.  We  took 
a i  to  be  the  typical  radiometer  noise  of  20  kR  at  all  data 
points  in  calculating  the  chi-squareds  in  Figure  23. 

Three  segments  from  the  data  base,  12,  3,  and  5,  appeare 
suitable  for  this  impulse -response  treatment.  Noise  and  the 
finite  intensity  resolution  of  the  data  limit  the  precision 
with  which  the  parameters  in  R  can  be  determined,  as  we  will 
find. 

RESULTS  --  EVALUATION 

The  cross-correlations  and  autocorrelations  of  I  are 
plotted  in  Figure  18a  and  b,  broken  out  into  groups  in  which 
the  aurora  was  fully  sunlit,  partially  sunlit,  and  not  sunlit. 
Tabulated  values  on  the  figure  give  t  at  the  peak  of  each 
curve  and  the  halfwidth  of  that  section  lying  above  the  normal 
ization  to  C(0)  at  1.00  (dashed  horizontal  line).  C(x) 
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Figure  18a.  Cross-correlation  functions  of  the 
15  aircraft  data  streams.  The 
tabulated  values  are  halfwidths 
above  the  normalization  at  1.0  of 
the  maximum  in  C(x). 
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dropping  below  this  zero-offset  level  implies  that  further 
relative  shifting  of  the  data  streams  will  result  in  a 
poorer  match  of  the  two  data  sets.  As  expected,  the  A's 
all  show  a  maximum  at  zero  lag;  no  subsidiary  peak,  which 
would  indicate  a  <10-sec  periodicity  in  the  aurora,  is 
evident  (except  perhaps  in  segment  15b). 

The  location  of  the  peak  in  the  cross-correlation 
function  should  not  depend  on  the  width  of  the  autocorrela¬ 
tion  function,  which  as  expected  increases  with  the  length 
of  the  data  segment.  We  interpret  the  location  of  this 
maximum  as  the  average  lag  between  SWIR  output  and  particle 
energy  input. 

(Note  that  R's  can  be  estimated  by  solving  the  above 
integral  equation  C  =  /RAdt^',  either  by  direct  fitting 
of  test  response  functions  or  by  the  usual  Fourier  transform- 
inverse  transform  procedure  applicable  to  convolutions; 
see  Ref  28.) 

These  peaks  are  compared  in  Figure  19a  with  the  lags 
estimated  previously  by  visual  inspection  of  the  data  streams. 
The  scatter  is  significant,  which  is  not  surprising  consider¬ 
ing  the  imprecision  of  the  visual  method.  (The  individual 
average  delays  from  the  cross  correlation  method  are  generally 
less  than  the  offsets  of  radiation  peaks  and  times  for  SWIR 
to  drop  to  baseline.)  Figure  19b  shows  a  weak  correlation 
between  the  lag  at  maximum  C  and  peak  of  the  energy  deposition 
altitude  profile.  Although  the  slope  of  the  straight-line 
equally  weighted  least-squares  fit  is  uncertain  and  the 
correlation  coefficient  is  low,  the  upward  trend  is  in  the 
direction  predicted  from  NO  chemistry.  Figure  19c  shows  that 
there  is  no  statistically  significant  dependence  of  the  ampli¬ 
tude  of  C(t)/C(0)-1  on  this  altitude.  The  cross  plot  in 
Fig  20,  in  which  the  lags  are  again  identified  by  segment 
number,  shows  no  significant  correlation  between  lag  at 
maximum  C  and  solar  illumination. 
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Figure  20.  Cross  plot  of  time  of  maximum  C(x 
illumination  of  the  auroral  profi 
FWHM  values  tabulated  in  Fig  18a. 


We  turn  now  to  deriving  the  parameters  of  the  model 
impulse  response  function  R  from  the  three  flight  segments 
in  which  aurora  turns  on  and  off  most  sharply.  Consecutive 
all-sky  views  of  segment  12's  sunlit  aurora  (the  perceived 
best  case)  are  in  Figure  21 ,  and  the  radiance  traces  are  in 
Figure  22.  The  center  of  the  small  open  circle  on  each  frame 
is  the  pointing  direction  of  the  -1/3°-fields  radiometer 
and  photometer.  As  the  aircraft  was  heading  north,  the 
instruments  pointed  26°N  of  the  magnetic  field  direction. 

Up  to  1125:26  a  relatively  low-level  (<25  kR  3914  A) 
diffuse  glow  lies  in  the  instrument  fields,  after  which  a 
bright  (>100  kR)  section  of  an  arc  that  had  formed  in  the 
west  sweeps  rapidly  eastward  across  them. 

Two  measured  quantities  enter  into  the  calculation  of  a 
a  set  of  Jp's,  K  (the  ratio  of  total  number  of  SWIR  to  3914- A 
photons  during  the  data  interval)  and  the  OH  radiance  (the 
hydroxyl  fundamental-band  background  level  within  the  radio¬ 
meter’s  bandpass).  In  practice,  the  goodness  of  fit  of  the 
computed  Jp's  is  sensitive  to  the  choice  of  both  these  para¬ 
meters.  The  bottom  trace  ir>  Fig  22  indicates  that  the  hydroxyl 
overtone  radiance  remains  nearly  constant  at  50  kR  during  the 
period  of  enhancement;  since  the  ratio  of  Av  =  1  to  av  =  2 
in-band  OH  photons  is  2.5:1  (Ref  27),  the  SWIR  pulse  would  be 
expected  to  be  superposed  on  a  -125  kR  OH  background. 

However  the  quiescent  level  of  the  SWIR  signal  derived  from 
the  actual  traces  is  about  355  kR,  which  we  adopted  as  the 
baseline  for  this  segment  (the  fits  using  125  kR  are  very 
poor).  Similarly,  the  value  of  K  that  we  had  derived  earlier 
(in  Ref  1,  3.1  photons/  photon)  is  both  low  compared  to  that 
from  other  segments  and  results  in  less-close  fits;  hence  we 
calculated  sets  of  Jp’s  using  also  K  =  4.6,  which  is  the 
ratio  at  the  peak  of  these  traces.  These  discrepancies  in 
absolute  infrared  signal  strength  may  indicate  some  error 
in  calibration  of  the  radiome'ter  in  the  segment. 
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Figure  22.  Radiometer  and  photometer  traces  of 
data  segment  12' s  aurora.  The 
supporting  data  at  top  and  bottom 
are  auroral  5577  A  radiance  (photo¬ 
meter  aligned  with  the  3914  A  photo¬ 
meter)  and  radiance  in  a  section  of 
the  hydroxyl  overtone-bands  airglow 
(wider-angle  radiometer  pointed  to 
zenith ) . 


Figure  23.  Chi-squared  resulting 
values  of  the  equilib 
segment  12. 
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Figure  24.  Pest-fits  of  SWIR  predicted  from 
model  excitation  function  to 
radiance  data,  segment  12. 
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Figure  23  plots  vs  xQ  for  n  between  0  and  2  for 
both  K's.  is  relatively  insensitive  to  n  and  xQ ,  even  for 

these  rapidly-changing  radiances,  as  neither  set  of  curves 
shows  a  well  defined  minimum.  When  n  =  0  (absolute  minimum  in 
X 2)  the  best  fit  is  at  xQ  =  1.4  sec  selecting  K  =.4.6,  and 
at  x 0  =  0.5  sec  for  K  =  3.1.  Jp(t)  in  this  case  reduces  to 
Kxq-1  /00t  I(t')  exp  f-(t-t")/x0l  dt'  (as  in  BRIM),  in  which 
the  average  lag  can  be  readily  shown  to  be  x0.  Figure  24 
plots  measured  and  predicted  SWIR  column  intensities  and  their 
deviation.  As  can  be  seen,  Jp(K  =  4.6)  --  that  is,  with  the 
more  plausible  energy  conversion  yield  --  gives  a  much  better 
fit;  its  standard  deviation  is  12  kR  (x  x  20  kR  noise), 
compared  to  39  kR  for  Jp(K  =  3.1).  The  1.4  sec  exponential 
time  constant  is  larger  than  the  0.8  sec  from  Fig  18a. 

In  segment  5  an  80  kR  3914  A,  half-sunlit  arc  some 
30  km  in  width  traverses  the  instrument  fields  closely  anti¬ 
parallel  to  the  aircraft's  motion  in  about  15  sec  (Figure  25). 
Figure  26  shows  results  from  the  calculation  of  Jp,  with  two 
n's  because  the  x^'s  plots  did  not  show  a  clear  difference 
between  minimums  at  the  best-fit  n,  xQ  pairs  (not  shown). 

The  closest  fit  is  reached  with  n  =  4.5,  xD  =  0.2  sec.  With 
n  forced  to  zero  xQ  is  2.5  sec  with  a  somewhat  poorer  fit; 
the  standard  deviation  is  41  kR  compared  to  26  kR  for  4.5, 

0.2  sec.  Clearly,  this  data  set  does  not  provide  good  resolu¬ 
tion  of  the  effects  of  predosing.  The  x0  derived  with  n  =  0 
is  again  about  twice  that  inferred  from  the  cross  correlation 
method. 

In  segment  5  the  aircraft  traversed  a  narrow,  weakening 
form  (Figure  27)  that  was  about  half  sunlit,  with  the  instru¬ 
ments  pointing  at  magnetic  zenith.  As  the  duration  of  SWIR 
enhancement  was  relatively  long  compared  to  the  expected 
x q  (15-20  sec),  it  represents  a  less  than  ideal  case 


Figure  25.  All-sky  photographs  of  data  segment  3's  aurora 


Figure  27.  All-sky  photographs  of  data  segment 
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Figure  28.  Fit  of  model  SWIR  radiance  function 
to  radiance  data,  segment  5. 


for  impulse-response  analysis.  The  best-fit  values  (Figure 
28)  are  n  =  0,  tq  =  2.1  sec;  again  td  is  about  twice  that 
derived  from  Fig  18a,  and  no  grow-in  factor  appears.  Note 
however,  that  the  fit  of  Jp  to  the  measured  J  is  poor  before 
and  after  passage  through  the  arc,  which  may  indicate  a  cali¬ 
bration  inaccuracy. 

CONCLUSIONS 

The  mean  delays  between  2.8  -  3.1pm-band  output  and 
particle  energy  input  inferred  from  cross  correlating  the 
two  auroral  data  streams  cluster  near  2  sec  and  vary  between 
0  and  7  sec  with  an  indication  of  increasing  with  the  height 
of  the  peak  in  the  energy  deposition  profile  between  -110 
-  140  km  altitude  (Fig  19).  Fits  of  three  arc  data  segments 
using  a  simplified  input-output  model  resulted  in  character¬ 
istic  exponential  SWIR  decay  time  constants  about  twice  those 
from  the  cross  correlation  peaks  with  no  resolvable  trend  with 
altitude  (the  best-fit  values  were  1-1/2,  2-1/2  and  2  sec). 
These  two-parameter  fits  suggest  that  the  emission  of  nitric 
oxide  vibrational  radiation  from  the  E  region  is  a  maximum 
immediately  after  impulsive  energy  deposition  (at  the  10^  - 
10^  ion  pairs/cm3  sec  level),  that  is,  that  most  of  the 
n2d  precursor  atoms  are  produced  directly  or  in  very  fast 
reactions. 

These  lag  times  are  within  their  resolution  consistent 
with  those  predicted  by  current  DoD  aerochemistry  models. 

(The  higher  total  yields  of  SWIR  nitric  oxide  overtone 
radiation,  however,  are  not  explained.)  Some  further  air¬ 
craft  data  segments  are  available,  and  other  analysis 
techniques  may  provide  better  estimates  of  the  time  delay 
and  geometric  offset  of  the  SWIR  emission. 
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